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SUMMARY 


This  is  the  third  and  final  report  of  a  three -year  study  on  compliant 
robotic  structures.  Such  structures  are  designed  of  continuously  flexible 
elastomeric  tubes  that  either  extend  as  bellows,  or  bend  transversely,  or  twist 
about  the  longitudinal  axis  when  subjected  to  internal  pressure.  The  motion  of 
such  a  tube  under  controlled  internal  pressure  depends  on  its  directional 
stiffness,  achieved  through  the  orientation  of  both  the  wall  corrugations  and 
reinforcing  fibers.  Tube  elements  in  series,  or  end-to-end,  make  up  a  robotic 
limb.  A  single  limb  may  be  employed  as  a  manipulator  arm.  Several  limbs  acting 
in  parallel  may  be  used  as  a  compliant  gripping  device  that  conforms  to  the 
object  being  manipulated;  and  multiple  limbs  in  parallel  may  be  used  as  legs  for 


walking  machines. 

^  -In- the  present  report,._Ha—have^  concentrate'^  on  the  analysis,  the  design,  and 
the  experimental  evaluation  of  the  above  mentioned  elements.  In  these  studies, 

^  wa  hove  tempere<f  out^conceptual  designs  with  a  knowledge  of  muscle  structure  in 

r 

selected  animals,  lie  .describe  a  computer-aided  control  scheme  applied  to  a 
demonstration  project:  a  manipulator  arm  made  up  of  bending  elements  in  series. 


This  50  cm  arm  successfully  manages  smooth,  open- loop,  pick-and-place 


manipulations  with  cycle  times  as  small  as  4.2  seconds  and  with  placement 


reproducibility  within  5%  of  the  target  point. 
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Chapter  I 

STATIC  BEHA^OR  OF  HIGH  FLEXURE  MANIPULATOR  ARMS 
Joseph  M.  Snyder  and  James  F.  Wilson 

INTRODUCTION 

A  high  flexure  manipulator  arm  is  shown  in  Figure  1.  The 
arm,  which  may  carry  an  end  payload,  is  comprised  of  several 
small  tubular  elements  and  is  moved  to  various  positions  by 
adjusting  the  internal  pressure  of  these  elements.  A  single 
element  of  this  arm  is  shown  in  Figure  2.  The  element  is 
fabricated  from  a  polymeric  material,  has  a  thin  corrugated 
bellows  section  for  part  of  its  circumference,  and  a  flat 
reinforced  strip  along  one  side.  Rigid  caps  are  provided  at 
each  end.  As  shown  in  Figure  2,  the  neutral  axis  of  bending 
of  the  element  is  located  in  the  flat  strip,  eccentric  from 
the  center  of  pressure  on  the  end  caps.  When  the  element  is 
pressurized  it  behaves  as  the  beam  shown  in  Figure  3,  where 
the  pressure  load,  F^,  rotates  to  follow  the  tangent  to  the 
elastic  line  and  is  offset  by  distance  d^.  The  purpose  of 
the  present  study  is  to  present  a  mathematical  model  to 
predict  the  large  static  deformations  of  a  single  element 
when  it  is  pressurized  slowly  and  lifts  an  end  payload,  W^. 

The  element  is  modeled  as  an  end  loaded  elastic  cantilever 


beam  for  which  self-weight  is  negligible  compared  to  the 
applied  loads.  Since  the  range  of  deformation  of  the  beam  is 
too  large  for  the  application  of  linear  beam  theory,  the 
beam  is  analyzed  using  the  exact  curvature  relations.  The 
shape  of  the  elastic  curve  of  a  beam  derived  from  the  exact 
curvatures  is  called  the  elastica* 
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The  elastica  has  been  studied  for  over  200  years.  Euler 
(1744)  performed  the  first  systematic  study  of  the  elastica. 
An  account  of  this  early  work  is  given  by  Timoshenko  (1953) . 
Love  (1944)  presented  solutions  for  the  "inflexional" 
elastica,  loaded  with  equal  and  opposite  forces  at  each  end 
as  shown  in  Figure  4a;  and  the  "non-inf lexional"  elastica, 
with  no  inflection  points,  as  shown  in  Figure  4b.  Timoshenko 
(1961)  presented  an  elastica  solution  for  a  cantilever  beam 
loaded  with  a  single  longitudinal  end  load  and  Hummel  and 
Morton  (1927)  gave  a  solution  for  a  single  transverse  end 
force.  Lau  (1982)  presented  a  solution  for  the  general  end 
load  problem  shown  in  Figure  5.  All  of  these  solutions  were 
non-explicit  analytical  expressions  using  elliptic 
integrals.  Additional  references  on  the  elastica  may  be 
found  in  the  literature  survey  by  Schmidt  and  DaDeppo  (1971) 
The  present  paper  deals  with  numerical  solutions  for  the 
elastica  of  a  cantilever  beam  subjected  to  general  end 
loads,  where  the  beam  models  a  pressurized  tubular  element 
with  an  arbitrarily  oriented  dead  weight  payload.  The 
results,  presented  in  nondimensional  form,  can  be  utilized 
in  two  ways:  to  develop  a  model  for  the  static  behavior  of 
the  entire  arm  and  to  model  the  arm  stiffness  needed  for  a 
dynamic  analysis  of  arm  motion. 
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MATHEMATICAL  MODEL 

Consider  the  cantilever  beam  with  end  loads  P^,  Q^,  and 
as  shown  in  Figure  5.  The  origin  of  the  axes  is  at  the 

tip  of  the  beam  and  the  loads  P^  and  remain  parallel  to 
these  axes.  The  basic  assumptions  of  the  element  model  are: 

i)  the  beam  is  weightless  and  all  loads  are  concentrated 
at  its  tip 

ii)  the  beam  material  is  linear  elastic 

iii)  the  beam  is  inextensible 

iv)  the  effects  of  Poisson's  ratio  and  transverse  shearing 
deformations  on  the  bending  deflections  are  negligible 

v)  unloaded,  the  beam  is  straight  with  length  L 

vi)  loaded,  the  curvature  of  the  beam,  dS/ds^,  is 
nonpositive. 

Based  on  these  assumptions,  the  differential  equation  of 
the  elastica  is 


de  ^ 


■  Qd^d  “  ^0 


where  E  is  the  modulus  of  elasticity  and  I  is  the  second 
moment  of  area  of  the  beam  cross  section.  With  the  following 
nondimensional  parameters 


®  “  L 


y  =  t: 


biL-a: 


equation  (1)  is  transformed  to  the  following  nondimensional 
form 

II  =  -Py  -  Qx  -  M  (3) 

When  equation  (3)  is  differentiated  with  respect  to  s  and 
the  following  relations  are  used 


equation  (3)  becomes 


^  =-Psin6  -  Qcos8 


This  equation  is  integrated  as 


d0  j / d0  X 

di 


=  I  (-Psi 


PsinS  -  QCOS0)  d0 


where  a  is  the  slope  of  the  elastic  curve  at  the  tip  of  the 
beam.  The  result  is 


§  =  P(cos0  -  cosa)  +  Q(sina  -sin0) 


When  equation  (7)  is  solved  for  d8/ds,  which  is 
nonpositive  in  this  problem,  the  following  equation  is 
obtained 


ds  = 


j2P(cos8  -  cosa)  +  2Q(sina  -  sin0)  + 


dx  =  COS0 

ds 

(4a) 

• 

dy  =  sin0 

ds 

(4b) 

M 
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(9) 


Equation  (8)  is  stated  more  concisely  as 

-d8 


ds  = 


^a  +  b  cos  ( 0+00 ) 


where 


a  =  -  b  cos  (a  +  0q)  (10a) 

b  =  (10b) 

0p  =  arctan|^j  ,  where  -7t<0(j<7i  (10c) 

Equation  (9)  is  integrated  over  the  length  of  the  beam  to 
obtain 


1-s 


_ d0 _ 

<Ja  +  b  COS0 


(11) 


where  ij)=0+0g.  Equation  (11)  is  a  non-explicit , closed  form 
solution  for  the  rotation  of  the  beam,  0,  at  the  arc  length 
coordinate  s.  The  nondimensional  coordinates  at  s  are 

COS(<J)-0n)  d0 

■  . — -  (12a) 

^ja  +  b  COS0 

®o 

sin(<j»-0f,)  d0 

y  =  -  ..  .  0'..^  (12b) 

J  -Ja  +  b  COS0 

®o 

Equations  (11)  and  (12)  may  be  expressed  in  terms  of 
elliptic  integrals,  but  this  provides  no  advantage  in 
obtaining  numerical  results  for  this  problem.  The  solutions 
remain  non-explicit  and  must  be  solved  numerically.  Given 
the  nondimensional  loads  P,  Q,  and  M  and  the  arc  length 
coordinate  s,  compatible  solutions  (0,x,y;a)  for  the  shape 
of  the  elastica  may  be  calculated  from  equations  (10)  - 
(12) .  The  loading  is  expressed  next. 
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LOADING  CONDITIONS 


The  static  behavior  of  a  single  element  is  modeled  by  the 
beam  with  the  loads  and  orientations  shown  in  Figure  6. 
and  N(,  are  the  force  and  moment  due  to  the  weight  of  the 
payload  at  the  tip  of  the  element.  Gravity  is  assumed  to  act 
at  the  constant  angle  P  as  shown.  is  the  pressure  load  at 
the  tip,  acting  at  eccentricity  d^  from  the  neutral  axis, 
see  Figure  2.  The  direction  and  magnitude  of  remains 
unchanged  as  the  element  rotates,  while  F^  rotates  to  match 
the  slope  at  the  tip  of  the  element.  varies  with  both  a 
and  the  configuration  of  the  payload. 

Define  the  following  nondimensional  loading  parameters; 


F.L^ 

^  El 

(13a) 

w  =  -ir 

(13b) 

(13c) 

=  T 

(13d) 

The  nondimensional  loads  of  equations  (10)  are 

P  =-F  cosa  -  W  cos() 

(14a) 

Q  =  F  sina  +  W  sinp 

(14b) 

M  =  N  +  F  d 

(14c) 

which  are  used  with  equations  (11)  and  (12)  to  calculate  the 
elastica  of  the  element  model. 
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ALTERNATE  EQUILIBRIUM  STATES 


An  elastica  problem  may  possess  more  than  one  possible 
solution.  Figure  7  illustrates  two  of  the  possible 
equilibrium  states  for  an  elastica  loaded  with  an  end  force 
Fp  If  F,  is  applied  directly  to  the  unloaded  beam,  the 
configuration  of  Figure  7a  will  result.  However,  the 
configuration  of  Figure  7b  will  result  if  the  load  F2  is 
applied  to  the  beam  before  F,  is  applied  and  removed  after 
F,  is  applied.  This  example  illustrates  that  the  loading 
history  must  be  defined  for  an  elastica  solution  to  be 
unique.  The  elastica  solutions  that  follow  are  unique 
because  for  a  fixed  initial  payload  the  internal  pressure 
loading  is  monotone  increasing,  which  is  consistent  with  the 
curvature  assumptions  used  in  the  derivation  of  equations 
(11)  and  (12) . 
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NUMERICAL  RESULTS 


The  prcjblem  shown  in  Figure  8  represents  a  vertical 
element  lifting  a  point  mass  of  weight  W^-  The  point  mass  is 
located  at  the  tip  of  the  element,  on  the  neutral  axis  so 
that  N^=0.  This  problem  was  solved  numerically  in 
nondimensional  form  using  equations  (10)-<12)  and  (14).  To 
solve  the  problem,  given  d,  F,  and  W,  equation  (11)  is  first 
solved  for  a  (s=0)  using  bisection  and  secant  iteration. 
After  a  is  calculated,  the  equations  are  used  in  an  explicit 
manner  to  determine  the  planar  coordinates  at  a  sufficient 
number  of  points  on  the  neutral  axis  to  describe  its 
deformed  shape.  Numerical  quadrature  is  used  to  evaluate  all 
integrals  in  the  equations.  The  computer  code  used  for  these 
calculations  is  in  Appendix  A.  Figures  9-16  show  the 
calculated  results  for  a  range  of  nondimensional  parameters 
representative  of  practical  experimental  elements.  Each  of 
these  figures  shows,  for  given  values  of  d  and  W,  the 
elastica  of  the  element  when  subjected  to  different  values 
of  the  nondimensional  pressure  load,  F.  For  clarity,  these 
figures  employ  the  fixed  X-Y  coordinate  system  of  Figure  8, 
normalized  by  the  element  length  L. 

Figures  9  and  10  show  the  elastica  with  no  dead  weight 
load  (W=0) .  As  expected,  for  a  given  level  of  pressure 
loading,  F,  the  curvature  and  the  tip  rotation  increase  as 
the  eccentricity,  d,  of  the  pressure  load  increases.  Figures 
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Fig.  8  Coordinate  system  and  loading  conditions  for  the 
numerical  examples  of  Figures  9-19 
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11-16  show  that  shortly  after  the  tip  of  the  element  rotates 
past  ninety  degrees,  a  relatively  small  increment  of 
pressure  load  causes  a  "flip-over",  or  large  change  in  shape 
of  the  elastica.  This  behavior,  not  observed  when  W=0 
(Figures  9  and  10) ,  becomes  more  pronounced  as  W  increases. 
Flip-over  is  of  special  interest  because  it  can  cause 
significant  problems  in  the  control  of  the  manipulator  arm. 

Figures  11  and  12,  for  which  W=2.5,  present  a  more 
detailed  view  of  the  flip-over  phenomenom.  Figure  11a,  for 
which  d=0.1,  shows  the  shape  of  the  elastica  as  F  varies 

from  0  to  15.  A  large  change  in  shape  occurs  as  F  increases 

from  11.8  to  11.9.  Figure  lib  shows  the  shape  of  the 
elastica  as  F  increases  from  11.8  to  11.9  in  increments  of 
one  hundredth.  The  element  still  changes  shape  in  an  abrupt 
manner  as  F  is  increased  using  these  finer  load  increments, 
the  change  occurring  between  F=11.84  and  F=11.85.  This 
abrupt  flip-over  also  occurs  in  the  elasticas  presented  in 
Figures  13  through  16.  Figure  12a,  for  which  d=0.2,  shows 
the  shape  of  the  elastica  as  F  varies  from  0  to  10.  A  large 

change  in  shape  occurs  as  F  increases  from  8  to  9.  Figure 

12b  shows  the  shape  of  the  element  model  as  F  increases  from 
8  to  9  in  increments  of  one  tenth.  The  element  changes  shape 
in  a  regular  manner  as  F  is  increased  using  these  finer  load 
increments. The  behavior  illustrated  by  Figure  12b  represents 


f  f  • 


a  tiransitlon  tetogeen  tlie  zegular  behavior  exhibited  in 
Figures  9  and  ID,  ozbere  W=<0,  and  the  abrupt  flip-over  which 
occurs  in  the  lemaiming  figures. 


Table  1  presents  the  pressure  load  at  the  onset  of  flip- 
over,  F^,  the  tip  rotation  at  the  onset  of  flip-over, 
and  the  change  in  tip  rotation  due  to  flip-over,  Aa,  for  a 
range  of  values  of  W  amd  d.  is  reported  to  the  nearest 
hundredth  and  Aa  is  computed  by  subtracting  from  the  tip 
rotation  calculated  for  When  W=2.5  and  d=0.2,  F^.^ 

is  reported  to  an  integer  valiae  and  and  Aa  are  not 
reported  because  the  change  in  shape  is  not  abrupt  and  hence 
these  quantities  cannot  he  precisely  identified. The  data  of 
Table  1  is  presented  graphically  in  Figures  17,  18,  and  19. 
These  data  show  that  flip-over  becomes  more  pronounced  with 
increasing  payload,  W,  and  with  decreasing  eccentricity,  d. 
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CHANGE  IN  TIP  ROTAT 


CONCLUDING  REMARKS 


Flip-over  is  the  result  of  interaction  between  the  payload 
W  and  the  pressure  load  F.  The  payload  always  acts  in  the 
same  direction  while  the  pressure  load  is  a  follower  load 
that  rotates  with  the  tip  of  the  elastica.  Before  flip-over, 
an  initial  payload  tends  to  return  the  elastica  to  its 
undeformed  shape.  When  the  element  is  near  vertical  (a  near 
0)  both  F  and  W  impose  a  tensile  load  on  the  element, 
increasing  its  stiffness  and  causing  deformations  to  be 
relatively  small.  When  the  tip  rotation  exceeds  ninety 
degrees,  the  vertical  component  of  F  opposes  W.  This  reduces 
the  restraining  effect  of  W  on  element  deformation  so  that 
small  changes  in  the  pressure  load  lead  to  flip-over. 
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Chapter  II 


DESIGN  AND  EXPERIMENTAL  EVALUATION  OF 
NEW  LIMB  ELEMENTS 

Jihad  N.  Ghattas  and  James  F.  Wilson 

INTRODUCTION 

Elastic  limb  elements  are  used  as  structural  components 
of  flexible  robotic  manipulator  arms.  Two  kinds  of 
elements,  bending  and  torque  elements,  can  be  linked  end  to 
end  to  form  an  entire  manipulator  arm.  See  Figure  1,  Such 
elements  incorporate  flexible  bellows,  and  can  be 
pressurized  independently  to  effect  the  motion  of  the 
manipulator  arm.  Pressurization  of  bending  elements  causes 
them  to  bend,  enabling  the  arm  to  lift  its  payload  (Wilson, 
1984a) .  Pressurization  of  torque  elements  causes  them  to 
rotate,  providing  twisting  action  for  the  arm  (Wilson  and 
Orgill,  1986). 

Most  robotic  arms  in  use  today  consist  of  hinged  or 
pivoted  joints  connected  by  rigid  links.  Flexible 
manipulators  utilizing  pressurized  elastic  elements  have 
several  advantages  over  such  rigid-link  arms,  including: 


bending  elements 


I 

I 

payload 


Figure  1:  Robotic  manipulator  arm  comprised  of 
bending  and  torque  elements.  Bending  elements  lift 
payload;  torque  elements  provide  twisting,  or 
"wrist”  action.  Pressure  to  each  element  is  inde¬ 
pendently  controlled. 


1)  Increased  durability  in  harsh  environments  due  to 
the  impact  resistance  and  ruggedness  of  elastomeric 

materials. 

2)  Reduction  in  arm  self  weight  due  to  the  remote 
location  of  the  power  source  (air  pressure  supply) . 

3)  Increased  speed  due  to  the  reduction  in  self  weight. 

4)  Increased  lifting  capacity  relative  to  self  weight. 


Disadvantages  include  an  increased  complexity  in  precisely 
controlling  the  arm^s  position  with  pressurized  air. 

There  are  several  previous  studies  that  are  relevent  to 
the  present  work.  For  instance,  the  instability  of  bellows 
subjected  to  internal  pressure  was  studied  by  Haringx 
(1952).  Wilson  (1984b)  evaluated  mathematical  models  for 
computing  the  overall  stiffness  for  bellows  of  various 
shapes.  Wilson  and  Orgill  (1986)  analyzed  uniformly  stressed 
orthotropic  cylindrical  shells  using  linear  theory.  Such  a 
model  may  be  used  to  describe  the  overall  deformation 
characteristics  of  the  torque  element. 

In  this  chapter,  the  mechanical  behavior  of  bending  and 
torque  elements  are  experimentally  evaluated.  Then  the 
measured  relationships  among  deflection,  internal  pressure, 


and  external  loading  are  compared  to  linear  theory. 
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DESCRIPTION  OF  ELEMENT  CONFIGURATIONS 


Three  types  of  elements  were  evaluated.  These  were: 


Bending  elements  with  uniform  size  of  bellows  gaps. 
Both  elements  with  one  and  two  bellows  (the  single- 
and  double-sided  elements)  were  evaluated. 

Bending  elements  with  bellows  gaps  of  optimum  size. 
Two  sizes  were  evaluated:  small  and  large. 

Torque  element. 


In  addition,  a  bending  element  with  six  bellows  (the 
"satellite"  configuration)  was  considered,  although 
quantitative  tests  were  not  performed  on  it  because  of 
irreparable  pressure  leaks. 

All  elements  were  made  of  type  A150  black  rubber.  The 
modulus  of  elasticity  of  this  rubber  was  previously  measured 
in  the  range  of  740  to  2330  psi.,  corresponding  to  the 
uncycled  and  cycled  tensile  specimens,  respectively  (Wilson, 
1986) .  End  fittings  were  aluminum  in  the  case  of  the 
bending  elements,  and  brass  in  the  case  of  the  torque 
element.  End  fittings  were  designed  to  allow  the  controlled 
intake  and  exhaust  of  air  while  preventing  leakage  of  air 


from  the  element. 


Bending  Elements  with  Uniform  Size  Bellows  Gaps. 

The  "single-sided"  bending  element  consists  of  a  solid 
slab  attached  to  a  corrugated  bellows.  See  Figure  2.  The 
bellows  consists  of  cylindrical  shells  with  radii  r,  and  rj 
connected  by  annular  plates.  A  detail  of  the  bending 
element  is  shown  in  Figure  3.  The  single-sided  element 
shown  in  Figure  2  has  equal  interior  and  exterior  gaps  (a  =  b 
in  Figure  3) .  The  thickness  of  the  bellows  wall  is  t. 

The  bellows  portion  of  the  bending  element  can  be 
modelled  as  a  cylidrical  tube  with  effective  wall  thickness 
tj,  radius  R  equal  to  the  bellows  mean  radius,  or 
(r|+r2)/2,  and  an  effective  modulus  of  elasticity 
E'  based  on  the  material  modulus  E  and  the  bellows  geometry 
(Wilson,  1984b) .  See  Figure  4a.  The  entire  bending  element 
is  modeled  as  a  composite,  cantilever  beam  of  length  L  with 
transverse  end  load  F^  and  moment  and  Mp.  See 
Figure  4b.  The  effect  of  longitudinal  end  loading  is 
neglected.  From  elementary  beam  theory,  the  transverse  tip 
deflection  is 

6  =  3^^  +  (Mo-Mp)2^ 

where  E  is  the  material  modulus  and  is  the  overall 
second  area  moment  of  the  composite  beam  section  consisting 
of:  (1)  the  slab  with  modulus  E  and  (2)  the  tube  model  of 

the  bellows  with  modulus  equal  to  its  effective  modulus  E' . 
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Figure  4;  Composite  beam  model  of  bending  element 
comprised  of  slab  and  tube  model  of  bellows  (4a) . 
4b:  cantilevered  beam  with  loading. 
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(2) 


is  calculated  by  first  locating  the  neutral  axis,  or 

\y  +  A,y 

V  sc  .  _ r 

^na  +  Aj 

where  the  cross-sectional  areas  of  the  cylinder  and  slab  are 


and 


A^  =  2itRtg 


(3) 


A,  =  b,h. 


(4) 


The  effective  thickness  t^  of  the  cylindrical  tube  is 


t  =  ^t 


(5) 


The  second  moment  of  area  about  the  section  neutral  axis  is  then 


I.  * 


b.h 


s  s 
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+  Ajy-y_ 


(6) 


The  moment  Mp  due  to  pressurization  is  given  by  the 

the  net  force  of  the  air  pressure  at  the  tube  end  acting  at 

a  distance  from  the  neutral  axis.  The  net  force 


acts  at  the  tube  center  and  is  the  product  of  the  pressure  and 
the  effective  area  at  the  tube  end,  or  pnR^.  Thus, 

Mp  becomes 

Mp  -  p-R^iyp-yp.)  < 

Program  LINDEFL.FOR,  listed  in  Appendix  B,  calculates 
element  deflections  using  Equations  (l)-(7)  where  E'  is 
based  on  the  Haringx  bellows  model. 

The  "double-sided"  bending  element  consists  of  two 
single-sided  elements  with  their  slabs  bonded  to  each  other. 
See  Figure  5.  The  two  bellows  are  on  opposite  sides  of 
the  attached  slabs  and  can  be  pressurized  independently  of 
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Fiaure  5:  Double-sided  bending  element  (a). 
Element  bends  up  when  pressure  in  lower  bellows  > 
nressure  in  upper  bellows  (b) ;  pressure  in  lower 
Snf  SppJr  belllws  equal  (e) ,  P«ssure  in  upper 
bellows  >  pressure  in  lower  bellows  (d) . 


each  other.  The  advantages  of  the  double-sided 
configuration  are: 


1)  The  element  can  bend  under  pressure  to  both  sides  of 
the  element's  neutral  axis  (instead  of  just  to  one 
side,  as  in  the  case  of  the  single-sided  element) . 

2)  Motion  can  be  better  controlled,  since  motion  due  to 
the  action  of  one  bellows  can  be  opposed  or  even 
reversed  by  the  action  of  the  opposing  bellows. 

Bending  Elements  with  Optimum  Size  Bellows  Gaps. 

The  bending  elements  with  the  optimum  size  of  bellows 
gaps  have  configurations  similar  to  the  single-sided 
elements  above,  except  that  the  interior  and  exterior  gaps 
(a  and  b  in  Figure  3)  are  not  equal.  Dimensions  a  and  b,  as 
well  as  radii  r,  and  rj  and  wall  thickness  t,  were  selected 
so  as  to  maximize  the  performance  of  the  element  (Ghattas, 
1988) .  The  performance  is  characterized  by  the  maximum 
degree  to  which  the  element  can  bend  toward  the  bellows 
side  without  a  closing  of  its  gaps.  The  element  must  also 
meet  deflection  and  payload  capacity  requirements. 

The  dimensions  of  the  small  and  large  elements  are  shown 
in  Figures  6  and  7.  These  elements  exhibit  a  maximum 
bending  angle  as  described  above  of  69  and  108  degrees, 
respectively.  This  is  a  significant  improvement  over  the 


Bending 


.04C  fyp 


maximum  angle  for  their  counterpart  elements  in  which  a  =  b 
and  the  maximum  bending  angle  was  45  degrees. 

I 

Torque  Element. 

The  torque  element  is  a  cylindrical  tube  with  bellows 
corrugations  at  a  helix  angle  of  70  degrees.  The  dimensions 
of  the  torque  element  are  shown  in  Figure  8 . 

The  end  rotation  of  the  torque  element  is  given  by 
Wilson  and  Orgill  (1986)  as 

=  i(i+ezz)Vs,  (8 

where  L  is  the  length  of  the  element,  R  is  the  mean  radius 
of  the  element,  is  the  longitudinal  strain,  and 
y.  is  the  shear  strain.  In  the  latter  reference,  the 
longitudinal  and  shear  strains  are  computed  for  various 
element  and  loading  parameters. 

Satellite  Bellows  Conf iguraton. 

The  bending  element  with  the  satellite  bellows 
configuration  consists  of  six  bellows  arranged 
circumferencially  (as  "satellites”)  about  a  flexible  central 
core.  Such  a  configuration  is  shown  in  Figure  9.  As  shown 
below,  six  satellites  is  the  number  which  yields  the  highest 
moment  due  to  pressurization  for  a  given  arm  ouside  diameter 


Figure  8:  Torque  element 


Figure  9;  Bending  element  with  satellite  bellow! 
configuration  (two  such  elements  connected  end  t( 


The  advantages  of  the  satellite  configuration  are  the 
same  as  those  of  the  double-sided  element  over  the  single¬ 
sided  element.  In  addition,  motion  is  not  limited  to  the 
plane.  Out-of-plane  motion  is  possible  because  the 
circumferencial  arrangement  of  the  bellows  enables  the 
element  to  bend  around  not  just  one  neutral  axis,  but  any 
neutral  axis  in  the  cross  section.  For  example,  in  Figure 
10,  if  bellows  1,  2,  and  3  are  pressurized  equally,  the 
element  will  bend  about  neutral  axis  a.  However,  if  bellows 
2,  3,  and  4  are  pressurized  equally,  the  element  will  bend 
about  neutral  axis  b.  Bending  about  any  intermediate  axis 
can  be  achieved  by  unequally  pressurizing  the  bellows,  or  by 
pressurizing  some  other  combination  of  bellows. 

The  optimum  number  of  satellite  bellows  was  determined 
to  be  six.  The  optimum  number  m  is  the  number  that  yields 
the  highest  bending  moment  for  a  given  outside  diameter 
OD  of  the  entire  cluster  and  a  given  internal  pressure  p. 

The  moment  M  due  to  the  pressurization  of  all  bellows  to  one 
side  of  the  neutral  axis  is  given  by  the  product  of  the  net 
force  due  to  the  air  pressure  p  (pressure  times  the 
effective  area  at  the  end  of  the  bellows)  and  the  distances 
of  those  bellows  from  the  neutral  axis.  Referring  to 
Figures  11  and  12  for  the  distances  of  the  bellows  centroids 
from  the  neutral  axis,  the  expression  for  moment  M  becomes 


Hi 


Figure  10;  Variable  neutral  axis  of  satellite 
configuration.  (a):  l,  2,  3  pressurized  yields 
bending  about  neutral  axis  a.  (b) :  2,  3,  4 
pressurized  yields  bending  about  neutral  axis  b. 
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Figure  12:  Element  with  m  satellite  bellows; 
m  is  odd.  Neutral  axis  is  throuqh  Section  which 
yields  highest  moment. 
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M  =  pnr^  X  sinf^^i+otj  i 

1=0 

Note  that  the  cluster  OD  determines  the  bellows  outside 


radius,  i.e.  (from  Figure  11) 


2  1+- 


i(n/m)  ] 


Substituting  this  expression  into  that  for  moment  due 
to  internal  pressure  gives  the  result 


9sin  (TT/m) 


[f 


sin (n/m) 


3  INT(m/2-l) 

X  sin[^^a]  (11) 

1=0 


where 


^  ^2Tt/m)  (for  m  even) 

a  =  or  Y 

|<2n/m)  (for  m  odd)  J 


and  INT  is  the  integer  function  that  trucates  after  the 
decimal  point. 

To  find  the  highest  moment  for  a  given  pressure  p  and 
cluster  outside  diameter  OD,  M  was  computed  for  various 
satellite  numbers  m.  The  results  tabulated  in  Figure  13 
show  that  for  six  satellite  bellows,  the  moment  will  be 
highest.  For  m>6,  the  moment  decreases  as  the  satellites 
become  infinitely  small  but  their  distances  to  the  neutral 
axis  approach  OD/2. 


1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18 


Number  of  Bellowe  m 


ure  13;  moment  in  element  with  unit  pressure  and 
side  diameter  2  as  a  function  of  number  of  satellite 
lows . 
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The  limb  shown  in  Figure  9  consists  of  two  bending 
elements,  each  with  six  comercially  available  satellite 
bellows  of  dipped  latex  rubber  that  surround  a  central 
core  made  of  a  smooth  rubber  tube.  Successive  pressuriza¬ 
tion  of  the  satellite  elements  did  produce  bending  about 
different  neutral  axes.  However,  because  of  the  inconsis¬ 
tency  of  the  bellows  wall  thicknesses,  irreparable  air 
leaks  occurred  even  at  pressures  as  low  as  10  psi  and 
experiments  of  a  quantatative  nature  could  not  be  performed 

DESCRIPTION  OF  THE  EXPERIMENTS 

The  purpose  of  the  experiments  was  to  test  the 
behavior  of  the  elements  and  to  compare  their  behavior  to 
the  results  predicted  by  linear  theory.  The  five  types  of 
elements  tested  were: 

1)  Uniform  bellows  gap  bending  element. 

2)  Double-sided,  uniform  gap  bending  element. 

3)  Small,  optimum  bellows  gap  bending  element. 

4)  Large,  optimum  bellows  gap  bending  element. 

5)  Torque  element. 

Experiments  performed  on  the  bending  elements  were: 


1)  Measurements  of  tip  deflection  under  internal 
pressurization.  Internal  pressure  in  the  bellows 
was  incrementally  increased  and  the  resulting  tip 
deflection  measured.  No  load  was  applied. 

2)  Measurements  of  load  capacity.  Load  at  tip  was 
incrementally  increased  and  the  pressure  required  to 
maintain  the  tip  at  zero  deflection  was  measured. 

3)  These  same  measurements  were  made  on  the  double¬ 
sided  element  where  the  initial  pressure  p^  in 

each  side  was  chosen  as  4  and  then  8  psi. 

Experiments  performed  on  the  torque  element  were: 

1)  Measurments  of  end  rotation  under  internal 
pressurization.  Internal  pressure  was  incrementally 
increased  and  the  resulting  end  rotation  measured. 

No  torque  was  applied. 

2)  Measurements  of  torque  capacity.  Torque  applied  at 
tip  was  incrementally  increased  and  the  pressure 
required  to  maintain  zero  end  rotation  was 
measured . 

3)  Measurements  of  end  rotation  under  torque  applied  at 
end.  Applied  torcjue  was  incrementally  increased  and 
the  resulting  end  rotation  measured.  Internal 
pressure  was  maintained  at  0  psi. 


The  experimental  apparatus  used  in  the  bending  tests  is 
shown  in  Figure  14.  The  tip  deflection  of  each  cantilevered 
bending  elements  was  measured  with  a  depth  micrometer  with 
an  accuracy  of  ±0.0005  in.  Loads  were  applied  by  suspend¬ 
ing  weights  from  the  tip  each  element. 

The  experimental  apparatus  used  in  the  torsion  tests  is 
shown  in  Figure  15.  Pure  torque  was  applied  to  this  element 
through  the  system  of  weights  and  pulleys.  End  rotation  was 
determined  by  measuring  the  deflection  of  a  rod  attached  to 
the  end  of  the  torque  element. 

All  elements  were  pressurized  with  laboratory  air  which 
passed  through  a  regulator  that  controlled  pressure  up  to  30 
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END  VIEW 


Figure  14 ;  Testing  apparatus  for  bending  and 
torgue  element  static  response  tests. 
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Figure  15:  Method  of  applying  pure  torque  to 

element  and  measuring  end  rotation. 
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EXPERIMENTAL  RESULTS  AND  COMPARISONS  TO  THEORY 


The  test  results  for  the  static  response  of  bending  and 
torque  elements  are  graphed  in  Figures  16-26.  Also  shown  on 
the  graphs  are  the  best-fit  lines  for  the  data  in  the  linear 
range  of  response  of  the  elements.  Their  slopes  are 
tabulated  in  Tables  1  and  2.  As  can  be  seen  in  the  graphs, 
most  tests  exhibited  essentially  linear  response  over  a 
certain  range,  after  which  the  response  became  nonlinear. 

The  results  are  also  compared  to  predicted  behavior.  To 
predict  element  behavior,  bellows  stiffnessess  were  computed 
using  the  Haringx  model  (Wilson,  1984b) .  The  bellows 
stiffnesses  were  used  to  calculate  the  relationships  among 
pressure,  deflection,  and  end  load  for  the  bending  elements 
and  among  pressure,  rotation,  and  end  moment  for  the  torque 
element.  These  are  also  tabulated  in  Tables  1  and  2,  where 
they  are  compared  to  the  experimental  values. 

Program  LINDEFL.FOR  was  used  to  compute  the  theoretical 
values  for  the  bending  elements  and  is  listed  in  Appendix 
B.  Torque  element  behavior  was  calculated  using  results 
from  Wilson  and  Orgill  (1986) .  Values  were  computed  for 
both  the  lower  and  upper  bounds  of  the  range  of  the  material 
modulus  E.  For  A150B  black  rubber,  experiments  show: 

E  =  740  psi  if  uncycled  and  E  =  2330  psi  after  cycling. 


Bending  Element 
Deflect  ion/ Pressure: 

Figures  16-19  show  the  vertical  tip  deflection  for 
bending  elements  subjected  to  internal  pressure.  The 
bending  elements  initially  responded  linearly,  but  the 
response  became  nonlinear  at  higher  pressure  and  deflection. 
This  is  to  be  expected  since  linear  theory  is  valid  only  for 
small  deformations. 

The  slopes  of  the  linear  response  ranges  are  listed  in 
Table  1,  where  they  are  compared  to  computed  values.  From 
equation  (1),  we  can  see  that  when  there  is  no  end  load  Fp 
and  no  end  moment  Mp,  the  vertical  tip  deflection  6 
becomes 

6  =  M,  (13) 

where  Mp  is  the  moment  due  to  the  internal  pressure  p 
acting  over  the  effective  area  at  the  bellows  end  at  a 
distance  from  the  bellows  center  to  the  neutral  axis  of  the 
beam.  The  computed  values  of  5/p  based  on  equations  (2)- 
(7)  and  (13)  are  listed  in  Table  1. 

Figure  17  shows  the  response  of  the  double-sided  element 
for  0,  4,  and  8  psi.  internal  pressure  in  the  opposing 
bellows.  No  significant  increase  in  element  stiffness  was 
observed  with  increasing  opposing  pressure. 

Pressure/ Load: 


Figures  20-23  show  the  load  capacity  of  the  bending 
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Figure  17 :  Experimental  vertical  tip  deflection 
as  a  function  of  pressure  for  double-sided, 
uniform-gap  bending  element,  with  no  load  at  tip, 
and  with  pressure  in  opposing  bellows  at  0,  4,  and 
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Figure  18 :  Experimental  vertical  tip  deflection 
of  small,  optimum-gap  bending  element  as  a 
function  of  pressure,  with  no  load  at  tip. 
Pressurization  and  depressurization  of  element 
shown . 
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Figure  19  :  Experimental  vertical  tip  deflection 
of  large,  optimum-gap  bending  element  as  a 
function  of  pressure,  with  no  load  at  tip. 
Pressurization  and  depressurization  of  element 
shown . 
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Figure  21  ;  Load  capacity  for  double-sided, 
uniform-gap  bending  element;  experimental 
pressure  required  to  sustain  load  with  zero  tip 
deflection.  Pressure  in  opposing  bellows  at  0,  4, 
and  8  psi.  Loading  and  unloading  of  element 
shown . 
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Figure  22:  Load  capacity  for  small,  optimum-gap 
bending  element;  experimental  pressure  required 
to  sustain  load  with  zero  tip  deflection. 
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Figure  23:  Load  capacity  for  large,  optimum-gap 
bending  element;  experimental  pressure  required 
to  sustain  load  with  zero  tip  deflection. 


elements,  represented  by  the  pressure  required  to  sustain 
the  load  with  zero  tip  deflection.  These  figures  also  show 
the  limits  of  linearity  for  these  elements. 


The  slopes  of  the  linear  response  ranges  are  listed  in 
Table  1  and  are  compared  to  computed  values.  When 
deflection  6  and  end  moment  are  zero,  equation  (1) 
becomes 


(14) 


3EIj, 

The  computed  values  of  p/F^  based  on  equations  (2) -(7)  and 
(14)  are  listed  in  Table  1.  The  ratio  p/Fg  is  computed 
for  a  single  value  rather  than  a  range  of  values  since  this 
ratio  is  independent  of  material  modulus  E. 


Torque  Element 

Figures  24-26  show  the  results  of  the  torque  element 
experiments.  The  slopes  for  the  linear  range  of  behavior 
are  listed  in  Table  2,  along  with  calculated  values. 

To  predict  the  behavior  of  the  torque  element,  the 
ratio  of  the  material  modulus  to  bellows  effective  modulus 
E/E'  for  the  bellows  was  calculated  using  the  Haringx 
bellows  model.  The  program  used  was  HARINGX. FOR,  which  is 
listed  in  Appendix  B.  For  the  torque  element  tested: 


This  value  was  then  used  with  the  results  of  Wilson  and 
Orgill  (1986)  to  predict  values  of  shear  strain 
as  a  function  of  bellows  material  and  geometric  parameters 
and  pressure  p  and  applied  end  torque  T.  y^^  was  then 
used  in  equation  (8)  to  obtain  the  desired  values.  Since 
our  experiments  showed  that  <<  1,  this  quantity  was 
neglected  in  subsequent  calculations  for  the  theoretical 
behavior  of  the  torque  element. 

Rotat  ion/ Pressure: 

The  measured  end  rotation-internal  pressure  data  are 
shown  in  Fig.  24.  Note  the  nonlinear  response  at  the  higher 
pressure  levels.  In  the  linear  range,  the  slope  (J)/p  is 
2.88  deg/psi. 

The  theoretical  ratio  0/p  for  this  torque  element, 
for  E/E'  =6.08,  is  based  on  the  following  results  given 
in  Fig.  2  of  Wilson  and  Orgill  (1986): 

V  -  (16) 

where  E  is  the  material  modulus,  t  is  the  bellows  wall 
thickness,  and  R  is  the  mean  bellows  radius.  With  equations 
(8)  and  (16)  with  =  0,  the  result  is 

I  =  (3.06)^  (17) 

Table  2  lists  0/p  for  lower  and  upper  bounds  of  E. 


Figure  2d:  Experimental  end  rotation  of  torque 
element  as  a  function  of  pressure,  with  no  torque 
applied  at  end. 


Pressure/T  orque: 


The  measured  pressure-torque  data  for  zero  rotation  of 
the  torque  element  are  shown  in  Fig.  25.  The  slope  p/T  is 
2.36  in-lb/psi  and  is  approximately  constant  up  to  p  =  10  psi. 

The  theoretical  ratio  p/T  for  E/E'  =6.08,  is  based  on 
Fig.  9  of  Wilson  and  Orgill  (1986),  or 

T  =  (3.68)pR^  (18) 

Solving  for  T/p, 

I  =  (3.68)R^  (19) 

where  T/p  is  independent  of  E.  Here,  R  =  0.5  in.  which  gives 
p/T  =  2.17. 

Rot  at  ion/T  orque: 

The  measured  rotation-torque  data  are  shown  in  Fig.  26. 
Here,  the  response  is  nearly  linear  for  positive  rotations  up 
to  10  deg  (rotations  following  the  helix  direction) ;  and  up 
to  20  deg  for  rotations  in  the  opposite  direction. 

The  shear  strain-torque  relationship,  based  on  Fig.  4  of 
Wilson  and  Orgill  (1986),  is 

=  (0.83)-^  (20) 

Substituting  into  (8)  and  solving  for  <>/T,  where  =  0, 
the  result  is 

Table  2  lists  0/T  the  lower  and  upper  bounds  of  E. 


Pressure  p  (pa) 


Figure  25:  Torque  capacity  of  torque  element; 
experimental  pressure  required  to  sustain  torque 
with  zero  «nd  rotation. 
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Table  1 

BENDING  ELEMENTS 
Summary  of  Experimental  and 
Theoretical  Behavior 


yi 


Test 

Peri  rmed 

Experimental 

Behavior 

Theoretical 
Behavior 
(E  =  740  to 
2330psi) 

Units 

deflection/ 

pressure 

6/p  =  0.884 

1.389  to  .4411 

in/psi 

pressure/ 

load 

P/Fq  =  12.0 

9.042 

psi/lb 

deflection/ 

pressure 

6/p  =  0.326 

.6280  to  .1994 

in/psi 

pressure/ 

load 

p/Fp  =  8.07 

9.042 

psi/lb 

deflection/ 

pressure 

6/p  =  0.481 

.9963  to  .3164 

in/psi 

pressure/ 

load 

P/Fq  =10.1 

8.084 

psi/lb 
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deflection/  6/p  =  1.31  2.671  to  .8483  in/psi 

pressure 

pressure/  P/Fq  =  6.64  3.526  psi/lb 
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Table  2 


TORQUE  ELEMENT 
Suimnary  of  Experimental  and 
Theoretical  Behavior 


Element 

Test 

Performed 

Experimental 

Behavior 

Theoretical 
Behavior 
(E  =  740  to 
2330psi) 

Units 

torque 

element 

rotation/ 

pressure 

0/p  =2.88 

6.95  to  2.21 

deg/psi 

pressure/ 

torque 

p/T  =2.36 

2.17 

psi/in-lb 

rotation/ 

torque 

0/T  =6.50 

15.1  to  4.77 

deg/in-lb 
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CONCLUSIONS 

Our  measured  ratios  deflection/pressure  (for  bending 
elements)  and  rotation/pressure  (for  torque  elements)  fell 
within  predicted  limits.  As  Tables  1  and  2  show,  these 
experimental  ratios  always  fell  between  the  predicted  upper 
bound  based  on  E  =  2330  psi  for  the  cycled  rubber,  and  the 
lower  bound  based  on  E  =  740  psi  for  uncycled  rubber.  In 
general,  the  measured  ratios  were  closer  to  the  theoretical 
values  based  on  the  upper  bound  modulus. 
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Our  studies  of  the  ratios  pressure/load  and  pressure/ 
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torque  involved  no  element  deformations  and  therefore  were 
independent  of  the  elastic  moduli.  For  three  of  the  four 
bending  elements,  the  measured  ratios  pressure/load  were 
higher  than  those  predicted  by  theory.  For  the  torque 
element,  the  measured  pressure/torque  ratio  was  within 
nine  percent  of  that  predicted  by  theory. 

Our  mathematical  models  that  predict  the  mechanical 
behavior  of  bending  and  torque  elements  rely  on  linear 
theory  and  a  knowledge  of  the  following  quantities:  Young's 
modulus,  the  shear  modulus,  and  Poisson's  ratio  for  the 
material;  the  loading  including  the  internal  pressure;  and 
element  geometry.  Our  experiments  herein  serve  to  validate 
these  mathematical  models  and  show  that  they  may  be  used 
with  confidence  in  designing  bending  and  torque  elements. 
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Chapter  III 


CONTROL  OF  FLEXIBLE  MANIPULATOR  ARMS 
James  F.  Wilson,  Zhenhai  Chen,  and  Rhett  T.  George,  Jr. 

INTRODUCTION 

Flexible,  corrugated  tube  elements  made  of  black  rubber  or  poly¬ 
urethane  are  designed  to  bend  when  subjected  to  internal  pressure.  A 
typical  configuration  is  shown  in  Figure  1.  Analyses  and  experimental 
tal  investigations  of  such  elements  were  carried  out  by  Wilson  and  Li 
(1987),  who  emphasized  the  control  of  two  elements  in  parallel  acting 
as  a  gripper.  Multiple  elements  in  series  such  as  shown  in  Figure  2 
can  also  be  employed  as  manipulator  arms.  Previous  analyses  of  such 
arms  were  made  by  Mahajan  (1985)  and  Palaniappan  (1986)  who  showed 
their  feasibility.  The  purpose  of  the  present  investigation  is  to 
evaluate  these  arms  experimentally. 

The  success  of  such  an  arm  as  a  fast  acting,  smoothly  operating 
manipulator  depends  on  its  pressure  control  system.  This  chapter  des¬ 
cribes  the  overall  experimental  system,  two  methods  of  pressure  con¬ 
trol,  and  typical  experimental  results.  The  experimental  system  con¬ 
sists  of  a  pneumatic  control  module,  an  I/O  board,  a  microcomputer, 
appropriate  software,  and  the  corrugated  tube  arm  elements.  The  first 
of  the  two  methods  of  element  pressure  control  is  based  on  timed 
sequencing;  the  valve-on  condition  in  which  air  is  admitted  continu¬ 
ously  until  the  desired  element  pressure  level  is  reached;  and  the 
valve-off  condition.  The  second  method  of  element  pressure  control  is 
based  on  pulsing  the  pressure  of  the  admitted  air.  Both  control  methods 
are  open-loop.  The  experiments  involve  using  these  methods,  together 
with  computer  programs  that  specify  different  pressure-time  histories 
for  the  arm  elements,  and  then  observing  the  speed  and  smoothness  of 
arm  and  gripper  motion.  In  a  typical  experiment,  the  gripper  at  the 
end  of  the  arm  grasps  the  payload  from  a  lower  platform,  the  arm  lifts 
the  load,  places  it  on  a  higher  platform,  and  returns  to  its  original 
position  without  mishap.  Experiments  show  that  the  fastest  time  for  a 
50  cm  long  arm  to  complete  such  a  cycle  is  approximately  five  seconds. 


80 


MANIPULATOR  ARM  AND  GRIPPER  DESIGN 


The  first  generation  arm  and  gripper  system  was  built  with 
seven  commercially  available  polyurethane  bending  elements.  The  out¬ 
side  dimensions  of  these  component  elements  (three  sizes)  are  shown 
in  Figure  1;  and  the  manipulator  design  is  shown  in  Figure  2.  The 
arm,  which  resembles  somewhat  the  trunk  of  an  elephant  in  both  its 
outward  form  and  manipulative  function,  is  clamped  at  the  top  and 
hangs  vertically  when  unpressurized  and  at  rest.  The  largest  ele¬ 
ments,  A-1  and  A-2,  are  at  the  clamped  end  and  these  elements  sus¬ 
tain  the  highest  arm  moments  when  the  pressurized  arm  with  the  pay- 
load  at  the  gripper  is  displaced  from  the  vertical  equlibrium  posi¬ 
tion.  In  terms  of  the  Cartesian  coordinate  system  defined  in  Figure 
2,  arm  elements  A-1,  A-2,  B-1,  and  C-1  are  aligned  to  bend  in  the 
x,y-plane;  but  element  B-2  is  aligned  to  bend  in  the  y,z-plane.  The 
latter  element  is  effective  in  avoiding  collisions  of  the  arm  with 
obstacles  in  the  x,y-plane.  Each  arm  element  is  a  distinct  pressure 
cell  and  has  its  own  air  supply  line  connected  to  the  pneumatic  con¬ 
trol  module.  The  two  opposing  elements  of  the  gripper,  however,  have 
a  common  air  supply. 

Consider  the  mechanics  of  this  manipulator  system.  When  air 
pressure  is  increased  in  a  particular  element,  its  bellows  portion 
expands  and  the  element  bends  about  its  neutral  axis  located  within 
its  flat  side.  See  Figure  1.  As  demonstrated  by  Wilson  and  Li  (1987) 
the  amount  that  a  given  element  bends  is  proportional  to  both  its 
internal  pressure  and  its  external  end  loads.  Thus,  for  the  manipu¬ 
lator  arm  to  function  properly,  it  is  apparent  that  the  internal 
pressure  supplied  to  each  of  its  elements  must  be  of  sufficient 
magnitude  to  overcome  the  opposing  forces  of  the  payload,  the  self 
weight  of  the  arm  and  gripper,  and  the  elastic  resistance  of  the 
element  material.  Mahajan  (1985)  and  Palaniappan  (1986)  studied  the 
interrelationships  among  these  system  parameters  and  their  results 
were  used  in  the  present  study  to  estimate  the  element  pressure 
levels  required  for  practical  manipulations  of  the  experimental  arm. 

A  typical  pick-and-place  maneuver  cycle  was  chosen  to  evaluate 


arm  performance,  especially  to  determine  the  shortest  possible  cycle 
times  consistent  with  smoothness  of  arm  motion.  A  photograph  of  this 
typical  maneuver,  taken  with  stroboscopic  light,  is  shown  in  Figure 
3.  Here,  the  gripper  is  pressurized  and  grasps  an  object  from  the 
lower  platform;  the  arm  elements  are  selectively  pressurized  to  lift 
the  object  to  an  upper  platform  while  avoiding  a  collision  with  it; 
the  gripper  is  depressurized  and  the  ooject  is  released;  the  arm 
elements  are  repressurized  to  raise  the  arm  somewhat;  and  then  all 
elements  are  depressurized  as  the  manipulator  returns  to  its  initial 
vertical  equlibrium  position.  Of  course,  other  maneuver  cycle.s  are 
also  possible,  including  the  retrieving  of  the  object  from  the  upper 
platform.  What  follows  are  descriptions  of  the  experimental  system 
used  to  achieve  these  maneuvers  and  discussions  of  typical  experi¬ 
mental  results. 

CONTROL  MODULE  AND  I/O  BOARD 

The  control  module  contains  12  solenoid  valves,  five  needle 
valves,  1/4  in.  diameter  flexible  tubing  and  connectors  for  the 
fluid  circuit,  and  the  necessary  electronics  for  interfacing  the 
control  module  with  a  microcomputer.  Each  of  the  five  arm  bending 
elements  and  also  the  gripper  has  an  independent  fluid  circuit  with 
two  solenoid  valves:  one  for  air  intake  control  and  one  for  exhaust 
control.  The  needle  valve  in  each  of  the  arm  element  fluid  circuits 
is  used  to  adjust  both  the  intake  and  exhaust  flow  rate  of  the  air. 
Figures  4  and  5  show  the  fluid  circuit  and  the  electronic  circuit, 
respectively,  for  a  bending  element.  In  a  typical  arm  maneuver,  the 
intake  air  pressure  for  each  of  the  arm  elements  is  preset  to  a  regu 
lated  value  of  40  psi.  However,  the  intake  pressure  for  the  gripper 
is  preset  at  a  higher  pressure  such  as  70  psi  to  insure  a  sufficient 
gripping  force. 

The  I/O  board,  a  commercially  available  Techmar  Lab  Tender, 
was  installed  in  an  IBM  pc-XT  microcomputer.  This  board  has  the  capa 
bilities  of  performing  four  main  functions:  analog  to  digital  con¬ 
version,  digital  to  analog  conversion,  digital  input-output  control, 
and  timing/counting.  The  digital  to  analog  conversion  and  counting 


Figure  3  -  Typical  pick-and-place  arm  maneuver:  A,  pickup 
point;  B,C,D,E,  lifting  cycle;  F,  placement  on 


upper  platform. 
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Figure  ^  -  Fluid  circuit  for  control  of  a  single  elemen 
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Figure  5 


Electronic  circuit  for  D/ A  control  of  air  intake 


and  exhaust. 
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features  of  this  Lab  Tender  board  were  utilized  to  develop  system 
software  for  arm  and  gripper  motion  control.  Two  succes.ral  software 
control  schemes  were  developed  and  these  are  discussed  next. 

SOFTWARE  CONTROL  SCHEMES 

It  became  quite  apparent  in  the  early  stages  of  this  research 
project  that  a  computer  was  needed  to  control  the  pressurized  air 
to  drive  the  arm.  For  instance,  when  an  operator  regulated  the  flow 
of  air  to  the  arm  elements  by  manually  pushing  buttons  to  activate 
the  solenoid  valves,  smooth  arm  motion  could  not  be  achieved.  Even 
after  some  practice  at  hand  control,  an  operator  could  not  decrease 
to  any  significant  extent  the  unwanted  arm  oscillations  in  attempt¬ 
ing  maneuvers  such  as  depicted  in  Figure  3.  Thus,  computer  control¬ 
led  schemes  were  devised  with  a  goal  of  achieving  such  a  maneuver 
cycle  in  10  seconds  or  less  with  motion  as  smooth  as  possible. 

Early  experiments  on  the  system  components,  the  solenoid  valves 
and  the  individual  bending  elements,  showed  that  their  response  times 
were  sufficiently  small:  12  milliseconds  for  the  solenoids  and  no 
longer  than  one  second  for  the  largest  bending  element  (type  A)  to 
fill  with  air  and  bend  a  full  90  degrees.  A  0.3  sec  response  time  was 
measured  for  elements  B  and  C  (Wilson,  1986). 

With  a  knowledge  of  the  solenoid  valve  and  bending  element 
response  times,  two  open-loop  software  control  schemes  were  developed 
to  achieve  the  pick-and-place  maneuver  of  Figure  3.  The  first  control 
scheme  is  the  timed  sequencing  method.  For  example,  Figure  6  depicts 
the  time  sequences  during  which  the  intake  and  the  exhaust  solenoid 
valves  are  open  to  drive  the  pressure  elements.  The  second  control 
scheme  is  the  pulsing  method.  Figure  7  shows  the  general  form  of  this 
control  signal  for  an  intake  or  exhaust  solenoid  valve  in  which  the 
valve  is  open  during  the  pulse  time  interval  and  closed  during  the 
pulse  gap  interval.  The  desired  arm  motion  is  achieved  by  fixing  the 
pulse  time,  the  pulse  gap,  and  the  number  of  pulses  for  each  solenoid 
valve  in  the  system.  To  implement  these  two  control  schemes,  the  D/A 
converter  and  the  timer  of  the  Lab  Tender  were  programmed  as  discus¬ 
sed  in  the  following  two  sections. 
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Figure  6  -  Timed  sequencing  for  a  typical  arm  maneuver. 
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PROGRAMMING  THE  D/A  CONVERTER 

To  operate  the  D/A  seetion  of  the  Lab  Tender,  the  following 
sequence  of  commands  are  used. 

1 .Select  a  voltage  to  be  output .  The  voltage  is  in  offset 
binary  which  means  that  binary  0000  0000  equals  -5.00  V,  binary 
1000  0000  equals  0.00  V,  and  binary  1111  1111  equals  +4.96  V. 
Therefore,  each  increment  of  one  bit  increases  the  voltage  by 
0.04  V.  The  base  address  of  t^e  Lab  Tender  board  is  8l6  and  the 
address  of  the  D/A  data  is  BAGE+5  =  8 1 6  +  5= 82 1  •  Thus,  a  command  to 
send  +4.96  V  out  is  {  in  BASIC 

10  OUT  821 ,255 

2 .  Enable  the  change  1  for  the  prev iously  selected  voltage  to 
be  output  on .  The  D/A  converter  has  15  channels  divided  equally 
among  two  multiplexers.  Multiplexer  0  controls  channels  0-7  while 
multiplexer  1  controls  channels  8-15.  The  eight  bit  binary  code 
is  set  up  so  that  D0-D2  selects  either  channels  0-7  or  8-15,  D3 
either  enables  or  disables  multiplexer  0,  D4  enables  or  disables 
multiplexer  1,  and  D5-D7  have  no  function.  Figure  8  shows  these 
bit  assignments.  Note  that  0000  1001=9  enables  multiplexer  0, 
channel  1,  and  that  0001  0011=19  enables  multiplexer  1,  channel 
10.  The  multiplexer  control  address  is  located  at  816+4=820.  The 
previous  two  statements  are  pr-igrammed  as: 
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20  OUT  820,9 
30  OUT  820,19 
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3. Disable  the  multiplexers .  This  must  be  done  between 
channel  changes  and  D/A  data  changes  or  glitching  may  occur.  To 
disable  multiplexer  0,  D3  must  be  0,  and  to  disable  multiplexer 
1,  D4  must  be  0.  Therefore,  outputting  0000  0000=0  disables  both. 
This  command  is: 


10  OUT  820,0 


To  illustrate  the  command  sequence,  consider 
program  that  outputs  +^1.96  volts  on  channels  0  and 


outputs 

0 

0 

• 

1 

1  volts 

10 

OUT 

821,255 

20 

OUT 

820 , 8 

30 

OUT 

820 , 0 

40 

OUT 

820  ,  17 

50 

OUT 

820,0 

60 

OUT 

821,0 

70 

OUT 

820,9 

80 

OUT 

820,0 

90 

OUT 

820,  18 

100 

OUT 

'  820,0 

annels  1  and  10. 


D/A  +H.96  V 
ENABLE  CHANNEL  0 
DISABLE  MULTIPLEXER 
ENABLE  CHANNEL  9 
DISABLE  MULTIPLEXER 
D/A  -5.00  V 
ENABLE  CHANNEL  1 
DISABLE  MULTIPLEXER 
ENABLE  CHANNEL  10 
DISABLE  MULTIPLEXER 


the  following 
9  and  then 


The  output  channels  of  the  D/A  multiplexers  are  sample/hold 
amplifiers  that  drift  to  their  most  negative  voltages  if  they  are 
not  refreshed.  Figure  9  shows  the  result  of  enabling  a  channel 
with  +4.96  V,  disabling  it  in  the  very  next  line  of  the  program, 
and  allowing  it  to  drift  to  its  most  negative  voltage.  This 
channel  was  connected  to  a  control  module  valve  and  the  time  and 
voltage  at  which  the  valve  closed  is  shown  as  the  circled 
coordinate  in  Figure  9. 
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PROGRAMMING  THE  TIMER 


operation,  the  program  must  first  set  the  Master  Mode  register, 
which  consists  of  a  16  bit  command  code  with  bits  assigned  as 
shown  in  Figure  10.  Since  the  IBM  PC  operates  with  an  8  bit  bus, 
the  command  must  be  an  input  command  for  two  consecutive 
statements.  Of  concern  here  is  bit  MM15,  which  is  for  scalar 
control;  and  bit  MMI^t,  which  is  for  Data  Pointer  control.  This 
program  uses  Binary  Coded  Decimal  so  MM15  is  1.  The  Data  Pointer 
control  is  used  to  enable  a  sequencing  mechanism  when  reading 
into  or  writing  out  of  specific  locations  in  the  timer. 

In  the  present  application,  the  internal  oscillators' 
frequency  of  1  MHz  was  divided  by  10,000  and  counted  by  one  of 
the  five  counters.  Each  counter  has  a  Mode  register,  a  Load 
register,  and  a  Hold  register.  The  Mode  register  establishes 
whether  the  counter  counts  up  or  down,  repetitively  or  once.  Its 
other  functions  are  not  used.  The  Load  register  contains  the 
value  that  the  counter  will  count  to.  In  this  application,  the 
frequency  coming  in  is  1  MHz/ 1 0 , 000= 1 00  Hz.  Thus,  if  the  maneuver 
cycle  time  is  set  at  15  sec.,  the  Load  register  contains  1500 
counts.  The  Hold  register  is  used  to  save  the  value  in  the 
counter  without  interupting  the  counting  process.  This  value  is 
read  by  the  computer  and  compared  to  the  preset  time  that  a 
particular  event  is  scheduled  to  occur. 

The  sequencing  referred  to  above  deals  with  the  manner  in 
which  the  Data  Pointer  accesses  the  Mode,  Load,  and  Hold 
registers.  If  MMIM  is  a  1  in  the  Master  Mode  Register,  the  Data 
Pointer  will  increment  through  the  registers  as  shown  in  Figure 
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The  bit  assignment  for  the  Counter  Mode  Register  is  shown  in 
Figure  12.  In  the  present  application,  the  counter  mode  was  set 
so  that  the  Count  Source  Selection  was  F5,  which  divided  the 
internal  oscillator's  frequency  of  1  MHz  by  10,000.  The  counter 
reloads  from  the  load  register  each  time,  it  counts  once  until 
armed  again,  and  it  counts  down.  The  counting  is  done  in  the 
Binary  Coded  Decimal. 


Counief  1  Mode  Reg 

I 

Counief  1  Load  Reg 

I 

Counter  1  Hotd  Reg 

I 

Counter  2  Mode  Reg 

I 

Counter  2  Load  Reg 

I 

Counter  2  Hold  Reg 

I 


I 

Counter  5  HoW  Reg. 


ELEMENT  CYCLE 


F igure  1  1 


IncreT.ents  of  the  data  pointer  in  an  elenen 


c  V  c  1  e 


Count  Sout 
oxxxx  = 

rce  Setection  — 
Count  on  Risir^g  Edge 

1XXXX  * 

Count  on  Falling  Edge 

xoooo  - 

TCN  -  1 

xoool  = 

SRC  1 

XOOlO  = 

SRC  2 

X0011  = 

SRC  3 

X0100  = 

SRC  4 

X0101  > 

SRC  5 

X0110  = 

GATE  1 

xom  * 

GATE  2 

XIOOO  • 

GATE  3 

X1001  *• 

GATE  4 

X1010  *- 

GATE  5 

X1011  » 

FI 

XitOO  *= 

F2 

Count  Control 

OXXXX  =  Disable  Sbecia'  Gale 
1XXXX  =  Enable  Special  Gaie 
XOXXX  *  Rdoaa  ifom  Loao 
xtxxx  e  Reload  ifom  Load  o'  Hold 
xxoxx  =  Couni  Once 
XX1XX  s  Count  Repetitively 
XXXOX  *  Binary  Cojnl 
XXX1X  -  BCD  Count 
XXXXO  Count  Down 
XXXXt  =  Count  U. 


CM16  CM14  CM13  CMt2  CM1 1  CWtO  CM9  CMB  CM7  CM6  CMS  CM4  CMS 


Gating  Cor>trol 
000  =  No  Gating 
OOt  =  Active  High  Level  TCN  *1 
OtO  »  Active  High  Leve'  GATE  N-  1 
011  *  Active  High  Level  GATE  N -  t 

100  »  Active  High  Level  GATE  N 

101  *  Active  Low  Leve'  GATE  N 

110  »  Active  High  Edge  GATE  N 

111  s  Active  Low  Edge  GATE  N 


Output  Control  > 

000  s  inaci'/c,  OjICj!  Low 

001  «  Activ;  Hign  le'Tiina’  Count  Pu'sr 

OiO  s  TC  T  igc'ec 

C*.  1  =  iiiega' 

100  *  inactive  O-'.pjt  Hign  impedance 

101  *  Active  Low  Te'ni’Hd'  Count  Fu'se 
no  =  ii-egai 

111  »  li.eca 


Figure  12  -  Bit  assignments  for  the  counter  mode  resister 


EXPERIMENTAL  RESULTS  AND  DISCUSSION 


For  the  manipulator  arm  and  gripper  defined  in  Figures  1  and  2, 
a  series  of  experiments  were  performed  to  achieve  the  pick-and-place 
maneuver  shown  by  the  photograph  of  Figures  3.  The  coordinate  system 
for  this  work  space  is  defined  in  Figure  13.  The  manuevered  object 
was  a  wood  cube  with  an  edge  dimension  of  2.54  cm.,  initially  orient¬ 
ed  with  its  sides  parallel  to  the  coordinate  planes,  with  its  eg 
located  at  position  X  =  Y  =  Z  =  0.  After  the  gripper  grasped  the  cube 
from  this  location  on  the  lower  platform,  the  arm  deposited  the  cube 
at  the  target  point  on  the  upper  platform,  position  X  =  47  cm. ,  Y  = 

49  cm,  and  Z  =  5  cm. ,  and  the  arm  returned  to  its  initial  vertical 
position.  The  total  time  for  this  maneuver  is  defined  as  the  cycle 
time  T.  The  performance  of  the  manipulator  and  the  two  control 
methods  were  evaluated  on  the  following  bases:  qualitative  observa¬ 
tions  of  the  smoothness  of  arm  motion;  quantatative  data  on  the  ampli¬ 
tude  and  frequency  of  arm  motion  at  the  end  of  the  maneuver;  measures 
of  the  deviations  AX,  AY,  AZ  from  their  respective  target  points; 
and  the  change  in  rotation  A0  of  the  cube  about  the  Y-axis,  from 
before  to  after  placement. 

The  first  series  of  experiments  employed  the  timed  sequencing 
control  method  with  the  basic  pressure  patterns  of  Figure  6.  A 
general  description  and  a  printout  of  a  typical  computer  program  for 
this  method  is  given  in  Appendix  C.  It  was  necessary  to  perform 
several  experiments  to  determine  the  needle  valve  settings  that  mini¬ 
mized  unwanted  arm  oscillations.  Qualitative  observations  of  these 
experiments  are  suummarized  as  follows.  The  smoothest  arm  motion 
occurred  when  the  needle  valves  were  set  to  greatly  restrict  the 
flow  of  air  and  the  maneuver  cycle  times  T  in  such  cases  were  as 
high  as  25  seconds.  As  the  needle  valves  were  opened  to  allow  for 
higher  air  flow  rates,  the  arm  moved  faster,  and  values  of  T  as  low 
as  5  seconds  could  be  achieved.  However,  this  fast  action  could  be 
achieved  only  at  the  expense  of  an  increase  in  unwanted  arm  oscilla¬ 
tions  during  the  maneuver  cycle.  Although  the  pick-and-place  senario 
was  achieved  for  many  trials  where  T  ranged  from  5  to  25  seconds, 
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there  was  still  a  need  to  achieve  smoother  arm  motion,  even  at  the 
higher  cycle  times. 

A  more  versatile  control  method  was  sought,  one  without  the 
need  for  needle  valve  adjustments  to  diminish  unwanted  arm  oscilla¬ 
tions.  This  lead  to  a  second  set  of  experiments  based  on  the  pulsing 
control  method.  A  general  description  and  a  printout  of  a  typical 
computer  program  for  this  method  is  given  in  Appendix  D.  In  this 
method,  each  solenoid  valve  was  activated  at  the  proper  preset  time 
during  an  arm  maneuver,  analogous  to  the  plan  in  Figure  6.  The  method 
employed  three  software  timers,  one  for  the  pulse  time  T1 ,  one  for 
the  pulse  gap  time  T2,  and  the  third  for  the  delay  time,  or  the  time 
required  to  grip  the  object  at  the  pickup  point  and  to  release  the 
object  from  the  grips  at  the  target  point.  In  programming  a  maneuver, 
these  three  time  intervals  may  take  on  different  values  and  this  ver¬ 
satility  allowed  for  the  design  of  maneuvers  with  smoother  arm  mo¬ 
tions  over  the  5  to  25  second  cycle  time  than  were  possible  for  the 
timed  sequenced  control  method. 

Typical  experimental  results  showing  the  accuracy  of  object 
placement  for  the  pick-and-place  maneuver  are  summarized  in  Table  1. 
In  each  of  the  two  programs,  the  time  intervals  for  the  maneuver 
cycle,  the  pulse,  and  the  gap  were  fixed.  For  each  program,  the 
change  in  the  object  rotation  A6  and  the  deviations  AX  and  AZ  from 
the  object’s  target  values  were  measured  in  repeated  trials.  Table  1 
gives  the  maximum  measured  A  &  and  the  maximum  percent  deviation 
from  the  target  coordinates,  based  on  ten  trials  for  each  program. 
Also  listed  is  the  maximum  amplitude  of  free  oscillation  A  which 
occurs  at  the  end  of  the  maneuver,  when  the  arm  is  depressurized  and 
swings  about  its  vertical  equlibrium  position. 

The  experimental  results  for  the  two  methods  of  arm  motion 
control  are  compared  and  summarized  as  follows. 

1.  If  the  gap  time  interval  T2  is  chosen  as  zero  for  the 
pulsing  method,  as  for  Program  2  of  Table  1,  then  the  two 
control  methods  are  the  same. 

2.  Both  control  methods  may  be  used  successfully  for  pick-and- 
place  manipulation  cycle  times  T  as  small  as  4.2  seconds. 

3.  The  pulsing  method  is  superior  to  the  timed  sequence  method 
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PROGRAM  STATIC 
REAL*8  ALPHA, XA,YA 

REAL*8  WEAR , PEAR , DEAR , EETA, PI , EPS , THETA 
REAL*8  A, E , E , RSQ , DCADRE , AERR , RERR , ERROR , F , FX , FY 
REAL*8  TEMP , TEMP  1 , TEMP2 , ALPHAl , ALPHA2 
REAL*8  C,X,Y,XT,YT 

INTEGER  II,III,IJ,IK,IL,IER,IFLAG,NCI,NC2 
INTEGER  IFILE.NL 
EXTERNAL  F,FX,FY 

COMMON  /SETl/  WEAR, PEAR, DEAR, EETA 
COMMON  /SET2/  E,KSQ, THETA 
COMMON  /SET3/  A,E,PI,IFLAG 

IFILE-I 

IFCIFILE  .EQ.  I)REWIND(IO) 

PI=3. IAI59265358979D+00 

EPS=. lD-06 

EETA-.OOD+OO 

WEAR-7.5D+00 

NL=20 

DEAR=0.20D+00 
AERR=O.OOD+O0 
RERR=. lD-09 
DO  1000  11=120,130,1 
PEAR=DFLOAT(II)/l.D+01 
IF(II  .EQ.  15)PEAR=1A. AD+00 
IFdl  .EQ.  20)PBAR  =  lA.5D+00 
PRINT* 

PRINT* 

PRINT  501 , WEAR, DEAR, PEAR 

PRINT*.'  INTERMEDIATE  COMPUTED  DATA’ 


TEMP1=-1 .OOD+00 

DO  1001  IJ=1,60,1 
ALPHA-DFLOAT (I J) *PI /20 . D+00 
CALL  COMP (ALPHA) 

IFdFLAG  .EQ.  1)G0TQ  1003 

TEMP2-DCADRE (F , A, E , AERR , RERR , ERROR , lER) - 1 . D+00 


IF(TEMP2  .GE.  0.0D+00)GOTO  1002 
1001  CONTINUE 

PRINT*, 'POSITIVE  LIMIT  COULD  NOT  BE  FOUND, STEPl’ 
STOP 

1003  ALPHA1-ALPHA-PI/20.D+00 
ALPHA2-ALPHA 
DO  1005  IJ-1,25,1 
ALPHA- . 5D+00* (ALPHA1+ALPHA2) 

CALL  COMP (ALPHA) 

IFdFLAG  .EQ.  1)THEN 
ALPHA2-ALPHA 
GOTO  1005 
END  IF 

TEMP2-DCaDRE(F. A. B. AERR. RERR. ERROR . 1ER)-1 .D*00 
IF(TEMP2  .GE.  0.  OD-^00)  THEN 
GOTO  1002  im 
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ENDIF 

ALPHAl'ALPHA 
1005  CONTINUE 

PRINT’-,  'POSITIVE  LIMIT  COULD  NOT  BE  FOUND. STEP2' 

STOP 

"1002  IF(DABS(TEMP2)  .LT.  EPS)G0T0  2000 
ALPHAl*0.0D+00 
ALPHA2-ALPHA 
DO  lOOA  IK-1,180,1 

ALPHA- ALPHAl- (ALPHA2-ALPHa1)  (TEMPI/ (TEMP2-TEMP 1) ) 

CALL  COMP (ALPHA) 

C 

TEMP-DCADRE (F , A, B , AERR , RERR , ERROR , lER) - 1 . D+00 
C 

IF (DABS (TEMP)  .LT.  EPS) GOTO  2000 
IF (TEMP  .GT.  0.0D*00)THEN 
TEMP2-TEMP 
ALPHA2-ALPHA 
ENDIF 

IF (TEMP  .LT.  0.0D+00)THEN 
TEMPI -TEMP 
ALPHAl -ALPHA 
ENDIF 

1004  CONTINUE 

PRINT’-, 'ALPHA  COULD  NOT  BE  FOUND  IN  180  ITERATIONS' 

STOP 

C 

2000  XA-DCADRE (FX , A, B , AERR . RERR , ERROR , lER) 

YA-DCADRE (FY , a, B , AERR , RERR , ERROR , lER) 

TEMP-DCADRE ( F , A , B , AERR , RERR , ERROR , I ER ) 

C 

PRINT  503.XA,YA, ALPHA, TEMP 
C 

PRINT-' ,  '  KSQ-  '  , KSQ ,  '  THETA-  ’  , THETA 
PRINT-',  'A-'  ,A.  'B-’  ,B 
PRINT-',  'E-'  ,E 
PRINT-' 

DO  1010  IJ=0,NL, 1 

C-DFLOAT  (I  J)  -'ALPHA/DFLOAT  (NL) 

C-A-^C/2.D+00 

X-DCADRE (FX, C , B , AERR , RERR , ERROR , lER) 

Y-DCADRE (FY , C , B , AERR , RERR , ERROR , lER) 

XT-XA-X 

YT-YA-Y 

PRINT  500,IJ,C,X,Y.XT,YT 
IFdFILE  .EQ.  1)WRITE(10,FMT-502)XT,YT 
1010  CONTINUE 
1000  CONTINUE 
C 

500  FORMAT (6X, 12 , 4X, 'C- ' ,D13 .6, 2X, 'X- ' ,D13.6, 2X, ’ Y- ’ ,D13.6, lOX, 
S’XT-' .D13.6,2X. 'YT-' ,D13.6) 

501  FORMAT (AX, 'WEAR- ' ,D13.6 ,2X, 'DEAR- ' ,D13.6,8X, ' PBAR- ' ,D13.6) 

502  FORMAT (2D16. 8) 

503  FORMAT (6X, 'XA- ' ,D23. 16, 2X, ' YA- ' ,D23. 16. 2X, 'ALPHA- ' ,D23. 16, 
SlOX, 'L-' ,D23. 16) 

STOP 

END 


SUBROUTINE  COMP (ALPHA) 

REAL’''8  alpha,  WEAR  ,  PBAR  ,  DEAR  ,  BETA,  TESTA .  TESTB 
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REAL*8  CA,SA,CB,SB,A,B,C,D,E,THETA,KSQ,PI,EPS,TEST 
INTEGER  I FLAG 

COMMON  /SETl/  WBAR , PBAR , DBAR , BETA 
COMMON  /SET2/  E,KSQ, THETA 
COMMON  /SET3/  A,B,PI,IFLAG 
EPS-l.D-09 
CA-DCOS (ALPHA) 

SA-DS IN (ALPHA) 

CB-DCOS (BETA) 

SB-DSIN(BETA) 

A—  (PBAR''CA+WBAR’'CB)  *2 .  D+00 
B—  (PBAR*SA+WBAR*SB)  *2 .  D+00 

C— A*CA-B*SA+  (PBAR*DBAR)  **2 
D=DSQRT(A*A+B*B) 

E«2.D+00/DSQRT(D+C) 

KS0“D/(D+C)+D/(D+C) 

THETA='DATAN2(B,A) 

TEST=D*DCOS (THETa) -A 
IF(ABS(TEST)  .GT.  EPS)GOTO  300 
TEST=D*DSIN (THETA) -B 
IF(ABS(TEST)  .GT.  EPS) GOTO  300 
A=-.5D+00*THETA 
B=.5D+00*ALPHA+A 
IF(KSQ  .LT.  1.0D+00)GOTO  100 
IF(DABS(a)  .GE.  (PI*.5D+00))GOTO  200 
IF (DABS (B)  .GE.  (PI* . 5D+00) ) GOTO  200 
TESTA=KSQ*DSIN (a) **2 
TESTB=KSQ*DSIN(B)**2 
IF (TESTA  .GE.  1.00D+00)GOTO  200 
IF (TESTS  .GE.  1 . OOD+00) GOTO  200 
100  IFLAG=0  ' 

RETURN 
200  IFLAG=1 
RETURN 

300  PRINT*, 'TEST  FAILED  IN  SUBROUTINE  COMP'  . 

STOP 

END 


FUNCTION  F(X) 

REAL*8  F,X,E,KSQ, THETA 
COMMON  /SET2/  E,KSQ, THETA 
F-1 .D+00-KSQ*DSIN (X) **2 
F-E/DSQRT(F) 

RETURN 

END 

FUNCTION  FX(X) 

REAL*8  FX,X,E,KSQ, theta 
COMMON  /SET2/E,KSQ, theta 
FX»  1 .  I>+00-KSQ*DSIN  (X)  **2 
FX-E/DSQRT(FX) 

FX-FX*DCOS (X+X+THETa) 

RETURN 

END 

FUNCTION  FY(X) 

REAL*8  FY.X.E.KSQ, theta 
COMMON  /SET2/E.KSQ. theta 
FY=1.-KSQ*DSIN(X)**2 
FY-E/DS0RT(FY) 
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APPENDIX  B 


Computer  Codes  for  Bending  Elements: 


LINDEFL. FOR 
HARINGX.FOR 


I 

s 

s 


§ 

I 


! 


i!* 

i 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCL'CCCCCCCCCrXCCCCCCCCCCLCCCCCCCCCCC 


C  ProcjrciiT,  LINDEFL.FOR  C 

C  C 

C  This  proyriiOi  computes  the  rei  ati  orishi  ps  between  de-f  1  ec  1 1  on .  C 
C  end  li3c<d,  end  nrtoment  tor  single  and  dot.tbl e-«ii ded  bending  C 
C  elemer.ts»  using  simple  beam  theoi  v  and  Haring^  model  -for  C 

C  bellows  sti-fFness.  C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


REAL  Mo , L , V AL  <  200 ) 


111 

FORMAT  (ID 

100 

WRITE  ’Enter ; 

001 

WRITE  (♦:,’<A)’> 
WRITE  (*, ’ (A) ’ ) ’ 

1 ) 

I NPUT 

DIFFERENT  VALUES 

WRITE  (*,’<A)’)' 

2) 

KEEP 

SAME  VALUES’ 

WRITE  (»,’(A)’)’ 
READ  <  * , i 1 1)  I 

3) 

QUIT 

PROGRAM’ 

wr  i  te  (  * ,  ’  <  ’  ’  I  fiput  E 
r  ead  <  * , * )  E 

;  ’ 

>  ’  ) 

(I.EQ.D  THEN 
WRITE  <*,’<ADD’ 

INPUT 

Si  :  ' 

READ  ’  (F6.4>  ’ 

)S1 

WRI  VEI  (^.’  (A'/’/' 

INPUT 

S2s  ’ 

RE  AD  (  *  •  ’  <  F  *) .  4  •  ’ 

)  S2 

WRITE  <)♦, ’(A)’)’ 

1  NFUT 

L  :  ' 

READ  (i',’(F6.4>’ 

)L 

WRITE  (*.’(A>'>’ 

INPUT 

Dl:  ’ 

REfiD  (  ♦• ,  '  ( F  P .  4  )  ’ 

)  D 1 

WRITE  D,’(A)')’ 

INFUV 

D2s  ’ 

READ  (*,’-;F6.4)’ 

)  D2 

WRITE  (*,’(A)’)' 

INPUT 

D3 :  ■' 

READ  (*,’(F6.4>’ 

)  D3 

WRITE  (*,’<A>')' 

INF  UT 

D4 :  ' 

READ  (*,  ’  (FE.4) 

)  1)4 

WR  I T  E  (*■,’<  A  )’>  ’ 

INPUT 

ALPHA  0 

READ  <»,’<F&.4)’ 

) ALPHAO 

WRITE  (iT.MA)’)' 

INPUT 

t  ;  •’ 

READ  (#,’(F6.4)’ 

>  t 

WRITE  ()ii,'<A)’)’ 

INPUT 

R  :  ’ 

READ  ’ <F6.4) ’ 

)R 

WRITE  OH, ’(A)’)’ 

INPUT 

Ro:  ’ 

READ  (*,’(F6.4)’ 

)  Ro 

WRITE  <*,’<A)’)’ 

INPUT 

H  !  ’ 

READ  <*, ’  <Ffc.4) ’ 

)H 

WRITE  (*,’(A)’)’ 

INPUT 

b  ;  ’ 

READ  (*,’(F6.4)’ 

)  b 

ELSE IF  (I.EQ.2)  THEM 
GOTO  007 

ELSEIF  (I.EQ.3)  THEN 
GOTO  iOOO 

ELSE 

WRITE  <*,  •'  <A>  •' )  •’  TRY  AIBAIN' 
GOTO  100 
END  IF 
EiTi=^0 . 00 
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8S 


<.>U  •• 


ii; 


. 5 

F.'HO^  ( 1 «  O'^Hy  h‘ )  /  ( 1 . 0'**H '  R ) 

RLAM=R4:T  /  b/b 

Qu=  <3.  OiK  <  1 . 0-V*V>  >  *  *0.  25/  <RLAM*  <  i 
Qi-  <3.0*  <  1. 0-V*V)  )*»<,..  25/  <RLAM*  <  1 
RLhMu- <SJNH(2.  ♦Qu)-*-SIN<2.  »au)  )/  <2 
RLAMi:^<SINH<2,  *Qi  ) +SIM<2.  *Qi  >  >  /  <2 
A” < i . -RH0*RH0+2. ♦RHO»RHO* ALOG  <RHG) ) *  *2*RLAMu 
lEi&=RHO*  ( 1 .  -RH0*RH0+2.  *ALC]16  <RHO)  )  **2*RLAni 
C- ( 1 . -RHDCRHO) * ( < 1 . -RHO*RHO) **2- (2. *RHD*ALOG  <RHO) ) **2> 
WRITE  (*, ’ <A) ' > 


0+H/ R ) ) **0.5 
,  0— H/R ) ) »  *0. 5 
,  *Qu> / <C0BH(2. 
-  *Qi ) / <C0SH<2, 


*Qu)  +C03<2.  *0.-0  ) 
*Qi  )  +-CC>&  2 .  *0i  )  p 


WR I TE 
WR I  TE 
write 
WR  I  TE 
WRITE 
WR  I  TE 
WRITE 
WR  1 1 E 
READ 


<*, 

(*, 

<*., 

(*, 

<* 

<* 

<* 

<  * 

<». 


(A) 
<A) 
<A) 
<A) 

<  A ) 
<A)  ’  > 
•;a)  •’ ) 


DO  YOU  WANT  TO: 


<A> ’  )  ' 

n  1 )  1 4 

fHEN 


i ) 


KEEP 

CALCULATE  CHI' 


IF  <iH.Et!.i> 

.  o 

ELSE  IF  <i4.EQ.2)  I  HEN 
CHJ~i  .  +  ^2.  »  <t./R)  )  .. 
END  IF 

FLAG1==3.  *R/4.  /t*CHI*  <1  . 


( I •  +  (H/R) >  * (A/C+DB/C) 


+H  R )  »  *  2  *  ( J  .  -  V  *  V )  /  RL AM  /  •;  B  /  R  >  »  » : 


.  1  •  '"FjHU*  *2) 


FLAe2^= 

FLAGl* ( ( I 

-RH0**2) **2 

-  (2. 

WR  I  TE 

(*, 

'  (A) '  ) 

WRITE 

(  ♦  , 

’  (A) ’  ) 

WR  n  £ 

(  »  . 

•’  <A)  '  ) 

WR  1  TE 

(  * , 

'  (A)  ■’  ) 

WRITE 

( * , 

'  (A) '  ) 

WR  I  TE 

( +  . 

'  (A)  '  ) 

CHOOSE  ONE 

WR  I  TE 

V*, 

’  (A>  ■'  ) 

WRITE 

’  (A>  ’  > 

1  ) 

OUTPUT 

DELT 

WRITE 

<*, 

'  (A) ' ) 

) 

OUTPUT 

DELT 

WRITE 

(*, 

’  (A) '  ) 

3) 

OUTPUT 

Fo/F 

WR  I  TE 

<*, 

’  (A) '  ) 

4) 

OUTPUT 

DELT 

Mq-=0’ 

WRITE 

(♦, 

•'  (A)  •'  ) 

5) 

OUTPUT 

DELT 

wr  i  te 

<*, 

'  (a) '  ) 

6) 

(  i) 

write 

<  * , 

•'  <  a )  ’  ) 

7) 

M  II 

(2) 

wr  j  te 

(♦, 

■'  <  a )  •' ) 

8) 

It  It 

(3) 

wr  ;i  t  e 

<*, 

’  <  a )  ’  ) 

9) 

It  It 

(4) 

wr  i  te 

(*, 

’  (a)  ■'  ) 

••• 

10) 

II  II 

(5) 

READ  < 

* , J 11 )  II 

FORMAT 

<F9.4) 

i/Fo  FOR  Ho-0  AND  R=0’ 

^/F  FOR  Fo=C>  AND 
FOR  DELTA-=<.*  AND  Miri=0' 
i  AT  MIDF’OINT  FOR  DEILTn  (END)  —  0  AND 

h/P  for  a  GIVEN  Fo  AND  Mo' 
double  eided  eleiiient’ 


IF  <  1 1 .  EQ.  1 .  or .  i  :i  .  eq- fc)  THEN 
PP=0.0 
Mo-0. 0 
Fo- i . 0 

ELSE IF  ( 1 I . EQ. 2. or . i i . eq. 7)  THEN 
F"o~0.  O 
Mo=0. O 
PP= 1 . 0 
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ELSE  IF  (I  I  .EQ.3.0R.  1 1  .EQ.  4.or  .  i  :i  .fe?q.  8.or  .  1  i  .  eq,9^  THEN 
Mo=0. 0 
PP= 1 . O 

ELSE  IF  (lI.ER.S.or.ii.eq. 10)  THEN 
WRITE  (♦,'<A)’)’  INPUT  Fos  ' 

READ  <*,  ■'  <F7.3)  =•  )  Fo 
WRITE  INPUT  Mo:  ' 

READ  <*, ' <F7.4) ’ )  Mo 

WRITE  (♦,'<A)’)'  INPUT  STARTING  PRESSURE:  ’ 

READ  <», ' (13) ’ ) IPS 

WRITE  <*,’(A)’)’  INPUT  ENDING  PRESSURE:  ' 

READ  ()«,  '  (13)  ’  )  IPE 

WRITE  (*.’(A)’)'  INPUT  DELTA  PRESSURE:  ’ 

READ  (*,  MI3)  '  )  IDP 
END  IF 

IF  (lI.EQ.S)  THEN 
00  11.4  I PF- 1  PS,  1  PE,  IDF 
PF-REAL(IPP) 

VAl.  (.  1  F'P )  ^VAl.  UE.  ( 1 1 ,  T  ,  F'P ,  D 1 ,  Dl;? ,  D3 ,  D4 ,  ALF FIAU ,  Rd  ,  R .  L ,  F'L  AG2 ,  Em , 

:l  E ,  Fo ,  Mo ,  YEAR ,  Y c  ,  RMUi.  T ,  Ap ) 

U4  CONTINUE 
ELSE 

VALM— VALUE  <  1 3. ,  t ,  F'P ,  D 1  ,  D2 ,  i)3 , 1)4  ,  ALPHAO .  Flo ,  R ,  L  ,  FL.AGl' ,  Efn ,  E ,  F o ,  Mo , 

1  '/EAR ,  V c  ,  RM(JL. T ,  Ap  ,  x  I L k o ) 

END  IF 

WRITE  (*  .  MA>') 

WRITE  (♦,"<A:«'') 

WRITE  (♦,’':A)') 
wr  3  <:  e  (  i  M  a )  '  ) 
wr  .3  to  (  ♦  ,  ■’  <«>  ■'  ) 
write 

i-f  (ii.eq.l.  or .  1 1  .  eq.  2 .  or .  i  i  .  eq .  3.  or .  ii  •eq.4.  or  .  i  i  •  eq-  ti)  then 
WRITE  (*,  •’  (A)  ■' )  •■■  Sir.gle-Sided  Eiemenf' 
el  se 

WRITE  v*'.,MA)')’  Double-Sided  £1  erner-.t ' 
endi  + 

WRITE  (i!,MA)') 

WRITE  <*, 151>Di,D2,D3 

Ibl  FORMAT!’  D1-',F9.4,’  D2=',F9,4,’  D3=’,F9.4) 

WRITE  (JK,  152)D4,t,R 

152  FORMAT  (’  D4=’,F9.4,’  t  =’,F9.4,'  R  =SF9.4) 

WRITE  <♦, 153) Ro, ALPHAO, L 

153  FORMAT!’  Ro=’,F9.4,’  ALPHA0=’ , F9. 4, ’  L=’,F9.4) 

WRITE  !#, 154>S1,S2,B,H 

154  FORMAT!’  S.t  =  ’,F9.4,’  S2-’,F9.4,’  Ei=’,F9.4,’  H^r’,F9.4) 

EPRIME=E./FLAG2 

WRITE  ( * , 1 55) E , Em , EPR I ME 

155  FORMAT!’  E  «’,F9.2, ’  Em=’,F9.4,’  Epr i me= ’ , F9. 4 ) 

WRITE  !», 157) YBAR, Yc,RMULT, Ap 

157  FORMAT  !’  Ybar=’ ,  F9. 4,  ’  Y'p=’,F9.4,’  E.!Hbar=’ ,  F9. 4,  ’  Ap=’ 

1,F9.4) 

Elbase  =  £♦•;>(  Ibase 

write!*,’ (’’  E*Iba«e  =  ’’, F9. 4) ’) Elbase 
WRITE  !*,15£>)CHI 

156  FORMAT  (’  CHI  =  ’,F9.4) 
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g 

6 

iS 

S 


P 

s 

i 


5 

i 


WRITE 
WRITE 
WRITE 
WRITE 
IF 


soo 


501 


503 


504 


5< 


>7 

,‘6 


.177 


<#. ' (A> 
<*,  ’  (A) 
U,  '  (A) 
(*,  '  (A) 


/, 


Fo=’ ,F4.1, 
, F9 . 4 ) 


Mo= ' , F4 . 


lio= '  ,  F4 , 


<  1 1  .EQ.  l.or  .  i  i  .eq.  6)  THEr4 
WIRITE  <*;.500>  Fo,MQ,F'P,VALri 
FORMAT  <'  FOR  THE  GIVEN: Fcj=*,F4. 

•  F-’ ,F4.2, ’  DEFLECTION  IBs  ',F9.4) 

ELSEIF  <I  I.  EQ.  2.c.r  .  i  i  .eq.7)  THEN 
WRITE  <ir  .501)  Fd,Mo,PP,  VALM 
FORMAT  <'  FOR  THE  GIVEN:’,./ 

-  P=',F4.2,/,’  DEFLECTION  IB: 

ELSEIF  <II.EQ.3.Dr.ii.eq.8)  THEN 
WRITE  (*,503)  DLETA,Mo,PP, VALM 

FORMAT  (•’  FOR  THE  GIVEN:',./,’  DELTA^’  r4.2,/,’  Ma=',F4,2. 

'  P=='  .F4.2,./,  '  FORCE  IS;  ',F9.4> 

ELSEIF  ( 1 1 .EO. 4.or . i i .eq.9)  THEN 
WP I TE  (  * , . 04 )  DELT  A , Mo , PF , VALM 
FORMAT  <  ■'  FOR  THE  GIVEN:  ’  ,  /,  ’  DELTA  (END) ,F4.2,  ./ ,  ' 

/  ,  ’  P=--' '  ,  F4 . 2 ,  /  ,  '  THE  MI  DPO I  NT  DEFLECT  ION  IS;  '  ,  F9 . 6 ) 

ELSEIF  (lI.EO.S.or . j i .eq. 10)  THEN 
WRITE  <*,505)Fo,Md 

FORMAT  <’  FDR  THE  GIVEN:  '  ,  / .  ’  Fo^-  ',F7.4.'  Mo- 
DO  506  I PP^ I PS, IRE. IDF 
WRITE  (  *  ,  507 >  I FP ,  VAl.  <  I PP ) 

FORMAT  (13,'  ',F6.4> 

CONTINUE 
END  IF 

(A)  ’  ) 

(A)  '  ) 

177) 

DO  YOU  WANT  TO  CONTINUE?  (1  FOR  YES.  2  FOR 
(II)’)  ICON 
THEN 


Mc.r^’  ,F4. 


' , F8. 4) 


1 0(T0 


WRITE  (*, 

WRITE  (*, 

WRITE  (*, 
FORMAT  (' 

READ  (*, 

IF  (ICON.EQ.i) 
GOTO  001 

ELSE 

GOTO  1000 
END  IF 
END 


140) 


i 

e 


FUNCT I  ON  VALUE  ( 1 1 ,  t ,  PP ,  D1  ,  D2 ,  D3 ,  D4 ,  ALPHAO ,  Ro ,  R ,  L. ,  '=’LAG2 ,  Em ,  E  ,  Fo , 
IMo,  YEAR,  Yc  ,  RMULT,  Ap,  >i  Ibsse) 

REAL  IDAR,  Ir,  Iw,  Is,  Ib,Mp,L,Mc. 

IF  (ALPHAO.lt. 3. 14159)  THEN 

ALPHA='2.  *  (ACOS  (CDS  (ALPHAO/ 2.  )  /R*Rd)  ) 

ELSE 

hLPHA^ALPHAO 
END  IF 

Se=(Em/E) *S1 
Te=( 1/FLA62) *t 
Ab=ALPHA*.R*Te 

Yb=2.  *R./ALPHA*S1N  (ALPHA/2.  )-Ro*COS( ALPHA/2.  ) +D2 

Ib--=2.  *R**3*Te*  (ALPHA-t-SIN  (ALPHA)  -4.  /  ALPHA*  ( 1 .  -COS  (ALPHA)  )  ) 


& 
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»i»(MW!iKiw»(nuuowjo\M»(n»uouGu<Mauo«H»u<niCiHiu(n»u(^^ 


Ar=Dl*D3 

Yr=.5*D3 

Ir= ( 1 . /12. ) »D1*D3»*3 
Aw=  2. *D4*  <D2-D3) 

Yw=^.5*  <D2+D3) 

Iw= ( 1. /6. ) *D4*  <D2-D3) **3 

AB=Sei|tS2 

Ys=-(.5)»S2 

Is~<  1.  /12.  )  i(tSe*32H(S2*B2 
A= Ab  +  Ar  +■  Aw+-As 

YBAR=  <1 .  /A)  t  (Yb*Ab+Yr#Ar+Yw*Aw+Ys)l(A©) 

XXX>;=Ib4  Ir  +  Iw+Is+AbUi  <YBAR-Yb)  »*2+Arf  <YBAR-Yr  )  *;»2 
IBAF:-XXXX-*-Aw*  (YBAR-Yw)  itt!ti2+AB*  (YBAR-Ys)  ♦*2 

>i  Ita&e=Ib4-Ir  +  Iw4-Is  +AbfYb#i)(24'  Ar!t(Yr**24-  Awi(iVw)(‘*2+  As*VsY»2 
I F  ( ALPHA .  L.T .  3 .  i  4 1 59  >  THEN 

Al-.5*ALPHAif^Rj(iR-RifRi*iS.TN<ALPHA/2.  )  *COS  <  ALPHA/2.  > 

A2---  i  D2-D3 )  *  <  D 1  --2 .  »  D4 ) 

Vc2-.5*  <D2+D3> 

XXX-2. *ft*R#R/3. /Al * ( 1 . -COS (ALFHA/2. ) *  *2>  »SJN  <ALPHA/2. ) 

YCl  =  XXX-R)|;COS(ALFHA/2.  )  +02 

Ap=Ai+A2 

wi"  i  to  t  ♦  5  ♦  y  '  Ap~ ’  ^  Ap 

I  c -  <  1 .  y  Ap )  H;  i  Yc.  1  i  +YC2 *  A2 ) 

ELSE 

A2-  <  D2-D3 )  *  <;  D  i  -2 .  *L.4 ) 

A3  -0 . 5 •*  ALPHA ♦  :R »;R 

VC3--4 . 0*R/  3.  O/ALPHAYSIN  <ALPHA/2.  O)  -Ro^CDS  <  ALPHA/  2.  0  <  +D2 
A4-— 0.5*Roik  (D1-2.0».D4)  *C0S  ( ALPHA./2. 0.« 
yC4— —  <  1 .  /3.  )  (ALF'HA/2.  /  +D2 

VC2=. 5* (02+03? 

Ap -A2+A3+A4 

Vc- m'  1 .  / Ap )  ♦.  •;  •t'C2 * A2+  YC3*;A3+  YC4  »: A4  > 

END  IF 

RHULT^E* JBAR 

Mp- ( Yc— YBAR ) *PF ♦ Ap 

i F  ( i 1 . eq . 6 . or . i i . eq . 7 . or .  i  i  . eq .  8 .  or  . 

1  i i . oq. 9. or . i  i  . eq. 1  > I  )  then 

y.  le-f-f  —  >i  Ibreise  ♦  2. 
el  ee 

>;Ie-f+  =  Ibar 
endi  P 

IF  (lI.EQ. l.OR. iI.EQ.2.DR. lI.EQ.S.or. 

1  i i . eq. 6. or . i i . eq. 7. or . i i . eq. 10)  THEN 

VALLlE=Fo*L**3/3.  /E/x  I eP P  +  (Mo+Mp )  *L*L./2.  /E/>:IePP 
ELSEIF  <II.EQ.3.or.ii .eq.8)  THEN 
VALUE  -•  -  (Mo+Mp)  *3.0/2.0/L 
ELSEIF  (II.EQ.4.or.ii.eq.9)  THEN 
VALUE=Mp H!L »L / 32 . 0/E /x  I  eP  P 
END  IF 

END 
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v*>yvvvi  x>c>c><yi<?c<«i 


i 

i 

I 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCrXCCCCCCCCCCCC 
C  Program  HARINGX.FOR  C 

C  Computes  bellows  stit  +  ness  according  to  Haring>!  model.  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
c  V«0. J 

c  b  ~ . 1 1 62 

c  R  =  (  . 5/cos <70. /  +  .5  )/2. 

c  t  --  .  1 

c  h  =  . 08 

RHO- < 1 . 0-H/R) / < 1 . 0+H/R) 

RLAM-R»:T/b./b 

Qu=  (3.  Ot  ( 1 . 0-V#V)  )  iK*0.  25/  (RLAM*  ( 1 . 0+H/R)  )  *  *0.  5 
Qi  =  <3. O* ( 1 . 0-V»V) ) **0.25/ (RLAM*  < 1 . 0-H/R> ) * *0. 5 

RLAMu--  (SINH  '2.  ^Qu)  +SIN  (2.  *0u)  )  /  <2.  VQu)  /  (COSH  (2.  »Qu)  +CDS  <2.  iKQuJ  ) 

RLAHi  =  (SINK  (2.  JKQi  )  4-SIN  (2.  *C>i  )  )  /  <2.  *Qi  )  /  (COSH  <2.  #Qi  >  4-COS  <2.  *Qi  )  ) 

A=  ( J  .  -RHD*RH042.  *RHD*RH0*AL0G  (RHO)  )  **2*RLAr'1u 

8B— RHO*  (  \  .  — F\H0*RH04-2 .  *ALOo<RHG.>  )  *  *2*RL  ANi 

C-  •;  1  .  -RHD*RHO)  *  <  •:  1  .  -RHQJfRHO)  **2-  (2,  *RHO* AL05  ( RHO)  )  **2) 


WR 

ITE 

<  * ,  ’ 

< 

A)  ’ 

)  ’ 

DO  I'GU 

I  WANT  TOs ' 

WR 

ITE 

< 

A)  =■ 

)  ' 

1  ) 

KEEP  CHl=i' 

WP 

ITF. 

<*,  •' 

( 

A)  ' 

)  ^ 

2) 

CALCULATE  CHI 

RE 

hD 

( * ,  * 

) 

14 

IF 

<14 

.  EQ, 

1 

f  T 

HEN 

CH 

1  -- 1 . 

0 

EL 

SEIF 

(14 

E.G'. 

2 ) 

THEN 

CH.T=--1 . 4  <2.  t  (b/  R>  )  /  ( 1 . 4  (H/R)  )  *  ( A/C4EijB/  C) 

END  IF 

FL  AGl  -=--.5.  *R/4.  /t*CHI  *  <  1 . 4H/R.>  *.*2*  ( 1 .  -V*V)  /RLAM/  <E/R)  ♦.*3/  <  i .  -RHO*  *2) 
FLAG2=FLH6i* ( ( 1 . -RHa**2>  **2- <2, *RHO*ALOS (RHO) ) **2) 

WRITE  (  *  ,  •'  ( ■'  '  Ey  Epr  i  me  -  '  '  ,  \  ) 

WRITE  (*,*)  "flag  2 
END 


I 


I 

e 


auou(^w^xuo^MU(^/^x^)<n)(^XMX)(Iu>>^A>laoc)eouffi»oool'Hvx)n(l^^ 


Pressure-Deflection 


Single-Sided  Element 


Dl  = 

0 . 5000 

D2= 

O . 0500 

D3= 

0.0500 

D4= 

O. OOOO 

t  = 

0. 0500 

R  = 

0.4000 

Ro= 

0 . 5000 

ALPHAO= 

5. 3600 

L  = 

3. OOOO 

Sl  = 

0.0000 

S2= 

0.0000 

B= 

0.0500  H= 

0. lOOO 

E  = 

740.00 

Em= 

0.0000 

Epr ime 

=  12.6557 

Ybar=  0.0618 

E*Ibase  =  0. 

Yp= 

7928 

0.4710 

E*Ibar= 

0.7168  Ap= 

0.5407 

CHI  =  2.4573 


FOR  THE  GIVEN: 

Fo=0. OO 
Mo=0. OO 
P= 1 . OO 

DEFLECTION  IS:  1.3890 


DO  YOU  WANT  TO  CONTINUE?  (1  FOR  YES.  2  FOR  N0> 


Si nal e-Si ded 

Element 

Dl  = 

0 . 5000 

D2= 

0 .  C»500 

D3= 

O.  05C>0 

D4  = 

C» .  OOOO 

t  = 

0 . 0500 

R  = 

0. 4000 

Ro= 

0.5000 

ALPHAO= 

5 .  36C»0 

L  = 

3 . 00C>0 

Sl  = 

Cl .  OC>OCi 

S2= 

0 .  OC>OC> 

B= 

0 . 0500  H= 

0. 1000 

E  = 

2330. 00 

Em= 

0 . OOOO 

Epr  i  me 

=  39.8483 

Yb  ar 

=  0.0618  Yn= 

0. 4710 

E*lbar= 

2.2571  Ap= 

Cl.  54Ci7 

E-^Ibase  = 

2.4961 

CHI  =  2.4573 


FOR  THE  GIVEN: 

Fo=0. OO 
Mo=0.00 
P=1 . OO 

DEFLECTION  IS:  0.4411 


DO  YOU  WANT  TO  CONTINUE?  <1  FOR  YES.  2  FOR  NO) 


Pres sure- He flection 


r 


p.  Double-Sided  Element 

■r. 


Dl  = 

0. 5000 

D2= 

0.0500 

D3= 

0.0500 

D4= 

O. OOOO 

t  = 

0.0500 

R  = 

O. 4000 

Ro= 

O. 5000 

ALPHAO= 

5.3600 

L  = 

3.0000 

Sl  = 

0.0000 

S2= 

0 . OOOO 

B= 

0.0500  H= 

O, lOOO 

E  = 

740.00 

Em= 

0.0000 

Eprime 

=  12.6557 

Ybar= 

0.0618 

Yp= 

0.4710 

E*Ibar= 

0.7168  Ap= 

0.5407 

■>'  E*Ibase  =  0.7928 

CHI  =  2.4573 


V 


FOR  THE  GIVEN: 

Fo=0.00 
Mo=0.00 
P= 1 . OO 

^  DEFLECTION  IS:  0.6280 

C  DO  YOU  WANT  TO  CONTINUE?  <1  FOR  YES,  2  FOR  NO) 

M! 


o 


i 


I 

r. 


Double-Sided 

Element 

Dl  = 

0. 5000 

D2= 

0. 0500 

D3= 

0.0500 

D4= 

0. OOOO 

t  = 

0. 0500 

R  = 

0.4000 

Ro= 

0. 5000 

ALPHAO= 

5 . 3600 

L  = 

3.0000 

Sl  = 

O. OOOO 

S2= 

0. OOOO 

B= 

0.0500  H= 

0. 1000 

E  = 

2330.00 

Em= 

0. OOOO 

Eprime 

=  39.8483 

Ybar 

=  0.0618  Yp= 

0.4710 

E*Ibar= 

2.2571  Ap= 

0.5407 

E*Iba5e  = 

2.4961 

CHI 

—  A«=; 

.73 

FOR  THE  GIVEN: 

Fo=0.00 
Mo=0. OO 
P=1 .OO 

DEFLECTION  IS:  O. 1994 


DO  YOU  WANT  TO  CONTINUE?  <1  FOR  YES.  2  FOR  NO) 


r> 

A 


119 


CL  cn  uj 


^ressure-ne flection 
(small,  optimum  bend.) 


Sinqle-Sided  Element 


Dl  = 
D4= 
Ro= 

Sl  = 

E  = 
Ybar= 


0.5000 
0.0000 
0.5230 
O. OOOO 
740.00 
O. 1914 


D2= 
t  = 

ALPHAO= 

S2= 

Em= 

Yp= 


0.3750 
0.0460 
5. 1070 
0.0000 
O. OOOO 
0.7941 


0.3750 

0.3980 

3.7500 

0.0777 


O. 1250 


Epr i me= 
E*Ibar= 


8. 1662 
2. 1831 


0.5133 


E*Ibase  = 
CHI  =  2. 


7.2930 

5795 


FOR  THE  GIVEN: 
Fo=0. OO 
Mo=0 . OO 
P=1 . 00 

DEFLECTION  IS: 


0.9963 


DO  YOU  WANT  TO  CONTINUE' 


(1  FOR  YES.  2  FOR  NO) 


Sinole-Sided  Element 


Dl  = 
D4= 
Ro= 

Sl  = 

E  = 
Ybar= 


O . 5000 
0 . OOOO 
0.5230 
0.0000 
2330. OO 
0.1914 


D2= 
t  = 

ALPHAO= 

S2= 

Em= 

Yp= 


0.3750 
0.0460 
5. 1070 
0.0000 
0.0000 
0.7941 


O. 3750 
0.3980 
3. 7500 
0.0777 


Eprime=  25.7123 
E*Ibar=  6.8739 


O. 5133 


E* I base  =  22.9632 

CHI  =  2.5795 


FOR  THE  GIVEN: 
Fo=0. OO 
t1o=0.00 
P=1.00 

DEFLECTION  IS: 


0.3164 


DO  YOU  WANT  TO  CONTINUE? 


<1  FOR  YES.  2  FOR  NO) 


'Mm*  i 


Load-Pressure 

(small  ontimum  bend.  el. 


rj 

•’  Sinqle-Sided  Element 


Dl  = 
D4= 
Fto= 

Sl  = 

E  = 
Ybar  = 


0.5000  D2= 

0.0000  t  = 

0.5230  ALPHAO= 
0.0000  S2= 

1000.00  Em= 

0.1914  Yp= 

9.8555 
2.5795 


E«Ibase  = 
CHI  = 


0.3750 
0.0460 
5. 1070 
0.0000 
o.c>ooo 

0.7941 


0.3750 

0.3980 

3.7500 

0.0777 


Eprime= 

E*Ibar= 


1 1 . 0353 
2.9502 


O. 1250 


0.5133 


FOR  THE  GIVEN: 
DELTA=???? 
lio=0 .  OCt 
P= 1 . 00 

FORCE  is:  -0.1237 


DO  YOU  WANT  TO  CONTINUE?  <1  FOR  YES.  2  FOR  NO) 


,4 


Sinale-Sided  Element 


Dl  = 
D4= 
Ro= 

Sl  = 

E  = 
Ybar= 


O. 7500 
0.0000 
0.7730 
O, 0000 
740.00 
O. 2835 


D2= 
t  = 

ALPHAO= 

S2= 

Effl= 

Yp= 


E*Ibase  =  34.6454 

CHI  =  2.5075 


Pres  su  re- De  ■'^lection 
(laree,  optimuin  bend,  el 


0.5620 
0.0460 
5.6400 
0.0000 
0.0000 
1 . 2268 


0.5620 

0.6010 

6.4540 

0.1128 


Epr i me= 
E*Ibar= 


4.6649 

9.5151 


O. 1720 


1 . 2936 


FOR  THE  GIVEN: 
Fo=0.00 
Mo=0.00 
P= 1 . OO 

DEFLECTION  IS: 


:.6709 


DO  YOU  WANT  TO  CONTINUE?  <1  FOR  YES. 


FOR  NO) 


Sinole-Sided  Element 


<;>.  7500 

D2= 

0.5620 

D3 

0. 0000 

t  = 

0.0460 

R 

0. 7730 

ALFHAO= 

5.6400 

L 

0. 0000 

S2= 

0.0000 

B= 

'330 .  OO 

Em= 

0. OOOO 

Epr 

0. 2835 

Yp= 

1 . 2268 

E*Ib 

D4=  0.0000  t 

Ro=  0.7730  ALPH 

Sl=  0.0000  S2 

E  =  2330.00  Em 

Yb ar=  0.2835  Yp 
E^Ibase  =  109.0862 

CHI  =  2.5075 


D3=  0.5620 
R  =  0.6010 
L  =  6.4540 
B=  0.1128 


O. 1720 


1 . 2936 


FOR  THE  GIVEN: 

F  0=0 . OO 
Mo=0 . OO 
P= 1 . 00 

DEFLECTION  IS: 


0.848' 


DO  YOU  WANT  TO  CONTINUE?  <1  FOR  YES.  2  FOR  NO) 


Loac-Pressure 

(larae,  optimum  bend,  el.) 


Si nql e-Si ded  Element 


Dl  = 
D4= 
Ro= 

Sl  = 

E  = 
Ybar= 


O.  7500 
0.0000 
0.7730 
0.0000 
1000.00 
O . 2835 


D2= 
t  = 

ALPHAO= 

S2= 

Em= 

Yp= 


0.5620 
0.0460 
5.6400 
0.0000 
0.0000 
1 . 2268 


D3=  0.5620 

R  =  0.6010 

L  =  6.4540 

B=  0.1128 
Eprime=  6.3039 
E*Ibar=  12.8582 


O. 1720 


1.2936 


E«lbase  =  46.8181 

CHI  =  2.5075 


FOR  THE  given: 
DELTA=?7?7 
Mo=0 . OO 
P= 1 . 00 

FORCE  is:  -0.2836 


DO  YOU  WANT  TO  CONTINUE' 


<1  FOR  YES, 


FOR  NO) 


'I 

i 


PROGRAM  I:  TIMED  SEQUENCE 


LINE 


DESCRIPTION 


Turns  off  any  channels  that  may  have  been 
left  on. 


10-20 


30-60 


70-100 


Defines  a  function  to  convert  Binary  Coded 
Decimal  into  decimal. 

Sets  the  Master  Mode  Register  (see  figure 
8).  Line  50  contains  the  low  byte  and 
line  60  contains  the  high  byte. 

Sets  the  Counter  Mode  Register  for  Counter 
#1  (see  figure  9).  The  low  byte  isin 
line  90  and  the  high  byte  is  in  line  100. 


110-130 


1M0-150 

160-180 


200-570 


Loads  2500  into  the  Load  Register.  Notice 
that  since  the  Data  Pointer  Increment  is 
enabled,  there  is  no  need  to  address  this 
register.  Once  the  counter  was  addressed 
in  line  80,  it  automatically  incremented 
to  the  load  register.  This  feature  makes 
it  very  easy  to  load  data  into  all  the 
counters  if  necessary. 

Disables  Data  Pointer  Increment 

This  command  loads  the  2500  from  the  Load 
Register  into  the  counter  and  arms  it. 
Arming  a  counter  starts  the  counting 
process.  Line  180  immediately  branches 
the  program  to  the  subroutine  which  grasps 
the  object. 

Branches  the  program  to  the  subroutine  at 
line  900  which  sends  a  command  to  save  the 
the  value  of  the  counter  in  the  Hold 
Register  of  the  counter  without 
interupting  the  counting  process.  The 
low  byte  contains  1/100's  of  a  second 
while  the  high  byte  contains  the  seconds. 

Checks  to  see  if  the  value  of  the  counter 
is  greater  than  the  time  the  next  event 
should  occur  (remember,  the  timer  counts 
down).  If  it  is,  it  goes  back  and 


S  >*ii  as 


refreshes  the  sample/hold  amplifiers.  If 
the  counter  is  less  than  the  time  of  the 
next  event,  it  goes  ahead  and  performs  it 


900-990 


Subroutine  for  checking  time.  Time  is 
saved  from  hold  register,  data  poniter 
is  put  back  to  the  hold  register,  and 
the  time  is  converted  from  BCD. 


1000-15270 


Subroutines  for  controlling  D/A  channels. 


iiCVilKKiCf! 


1  GOSUB  15000 
10  REM 


15000  'TURN  ALL  VALVES  OFF 

DEFINE  BCD-DECIMAL  CONVERSION 
FNDCONV (X) =X-6* (X\16) 

SET  MASTER  MODE 


DEF  FNDCONV(X)=X-6* (X\16) 

REM  SET  MASTI 

OUT  825,23 
OUT  824,0 
OUT  824,128 
REM  SET  COUNTER  MODE 

OUT  825,1 
OUT  824,16 
)  OUT  824,15 

)  REM  SET  LOAD 

)  OUT  824,0 
)  OUT  824 , 37 

)  REM  DISABLE  5 

)  OUT  825,232 

)  REM  LOAD  AND 

)  OUT  825,97 
)  GOSUB  1000 
)  GOSUB  900 

)  IF  024.8  THEN  GOTO  180 
)  GOSUB  2000 
)  GOSUB  900 

]  IF  024.6  THEN  GOTO  210 
)  GOSUB  3000 
)  GOSUB  900 

)  IF  024.3  THEN  GOTO  240 
)  GOSUB  4000 
)  GOSUB  900 

)  IF  0  24  !  THEN  GOTO  270 
)  GOSUB  5000 
)  GOSUB  900 

)  IF  023.8  THEN  GOTO  300 
)  GOSUB  6000 
)  GOSUB  900 

)  IF  021!  THEN  GOTO  330 
]  GOSUB  7000 
')  GOSUB  900 

}  IF  019.1  THEN  GOTO  360 
)  GOSUB  9000 
3  GOSUB  900 

)  IF  018.9  THEN  GOTO  390 
)  GOSUB  8000 
)  GOSUB  900 

)  IF  018.5  THEN  GOTO  420 
3  GOSUB  10000 
3  GOSUB  900 

L  IF  018.3  THEN  GOTO  450 
I  GOSUB  10500 
i  GOSUB  900 

3  IF  017!  THEN  GOTO  46  2 
3  GOSUB  11000 
3  GOSUB  900 

3  IF  016.5  THEN  GOTO  480 
3  GOSUB  12000 


LOAD  FOR  25  SECONDS 


DISABLE  SEQUENCING 


BEGINNING  OF  SEQUENC 


520 

530 

540  ' 

550 

560 

570 

580 

900 

901 
910 
920 
930 
940 
950 
960 
965 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
2000 
2010 
2015 
2020 
2030 

2040 

2041 

2042 

2043 

2044 
2050 
3000 
3010 
3020 
30  30 
304  0 
3050 
3060 

3070 

3071 

3072 

3073 

3074 
3080 
4000 
4010 
4020 
4030 
4040 
4050 


GOTO 


GOSUB  900 
IF  016!  THEN  GO 
GOSUB  13000 
GOSUB  900 
IF  09!  THEN  GOTO 
GOSUB  15000 
END 
REM 
REM 

OUT  825,161 
REM 

OUT  825,17 
REM 

A=INP (824 ) 

B=INP (824 ) 
B=FNDCONV(B) 

A= . 01*FNDCONV ( Aj 
C=A+B 
RETURN 
REM 

OUT  820,0 
OUT  821,255 
OUT  820,8 
OUT  820,0 


SUB  FOR  CHECKING  TIME 
SAVE  COUNTER  1 


DATA  POINTER  HOLD 


READ  COUNTER 


OUT  82C 
OUT  821 
OUT  82C 
OUT  82C 
RETURN 
REM 

OUT  821 
OUT  82C 
OUT  82C 
OUT  82C 
OUT  82 C 
REM 

OUT  82] 
OUT  82 ( 
OUT  82( 
RETURN 
REM 

OUT  821 
OUT  82( 
OUT  8  2( 
OUT  82( 
OUT  82( 
OUT  82( 
OUT  82( 
REM 

OUT  82] 
OUT  82( 
OUT  82( 
RETURN 
REM 

OUT  82] 
OUT  82C 
OUT  82( 
OUT  82( 
OUT  82C 


SUB  FOR  TURNING  ON  CHANNEL  8 
TURN  MULTIPLEXERS  OFF 
4.96  VOLTS 
ENABLE  CHANNEL  8 
DISABLE  CHANNEL  8 


255 

8 

0 

,22 
,  0 

LINES 

,  0 


SUB  TO  TURN 
'5.00  VOLTS 
' ENABLE 
'DISABLE 


KEEP 


CHANNELS 


FOR  TURNING 


CHANNELS  8,22,18 


255 

8 

0 

22 

0 

18 

0 

LINES 

0 


'ENABLE  8 
' DISABLE 


KEEP 


TURN 


8,22,20,18 


V  *•"  “V  *-• 


4060  OUT  820,20 
4070  OUT  820,0 
4080  OUT  820,18 

4090  OUT  820,0 

4091  REM  LINES  TO  KEEP  16  OFF 

4092  OUT  821,0 

4093  OUT  820,16 

4094  OUT  820,0 
4100  RETURN 

5000  REM  SUB  TO  TUI 

5010  OUT  821,255 

5020  OUT  820,8 

5030  OUT  820,0 

5040  OUT  820,22 

5050  OUT  820,0 

5060  OUT  820,20 

5070  OUT  820,0 

5080  OUT  820,18 

5090  OUT  820,0 

5100  OUT  820,16 

5110  OUT  820,0 

5120  RETURN 

6000  REM  SUB  TO  TUI 

6010  OUT  821,255 

6020  OUT  820,8 

6030  OUT  820,0 

6040  OUT  820,22 

6050  OUT  820,0 

6060  OUT  820,20 

6070  OUT  820,0 

6080  OUT  820,18 

6090  OUT  820,0 

6091  OUT  820,16 

6092  OUT  820,0 

6093  OUT  820,10 
6100  OUT  820,0 
6110  RETURN 

7000  REM  SUB  TO  TUI 

7010  OUT  821,0 

7020  OUT  820,20 

7030  OUT  820,0 

7040  OUT  820,16 

7050  OUT  820,0 

7060  OUT  821,255 

7070  OUT  820,8 

7080  OUT  820,0 

7090  OUT  820,22 

7100  OUT  820,0 

7110  OUT  820,18 

7120  OUT  820,0 

7130  OUT  820,10 

7140  OUT  820,0 

7150  OUT  820,17 

7160  OUT  820,0 

7170  OUT  820,21 


SUB  TO  TURN  ON  8,22,20,18,16 


SUB  TO  TURN  ON  8,22,20,18,16,10 


SUB  TO  TURN  OFF  20,16  AND  TURN  ON  8,22,18,10,17,21 


7180 

OUT 

820,0 

7181 

REM 

LINES 

TO  KEEP 

16 

OFF 

7182 

OUT 

821,0 

7183 

OUT 

820,16 

7184 

OUT 

820,0 

7190 

RETURN 

8000 

REM 

SUB 

TO 

TURN  OFF  8  AND  TURN  ON  22,18,10,21,9 

8010 

OUT 

821,0 

8020 

OUT 

820 , 8 

8030 

OUT 

820 , 0 

8031 

REM 

LINES 

TO  KEEP 

16 

OFF 

8032 

OUT 

820,16 

8033 

OUT 

820 , 0 

8040 

OUT 

821,255 

8050 

OUT 

820,22 

8060 

OUT 

820,0 

8070 

OUT 

820 , 18 

8080 

OUT 

820  ,0 

8090 

OUT 

820 , 10 

8100 

OUT 

820 , 0 

8130 

OUT 

820 , 21 

8140 

OUT 

820,0 

8150 

OUT 

820 , 9 

8160 

OUT 

820 , 0 

8161 

REM 

LINES 

TO  kee: 

16 

OFF 

8162 

OUT 

821,0 

8163 

OUT 

820 , 16 

8164 

OUT 

820,0 

8170 

RETURN 

9000 

REM 

SUB 

TO 

TURN  OFF  17  AND  TURN  ON  22,18,21,8,10 

9010 

OUT 

821,0 

9020 

OUT 

820  ,  17 

9030 

OUT 

820,0 

9031 

REM 

LINES 

FOR  16 

9032 

OUT 

821,0 

9033 

OUT 

820,16 

9034 

OUT 

820,0 

9040 

OUT 

821,255 

9050 

OUT 

820 , 10 

9060 

OUT 

820 , 0 

9070 

OUT 

820,22 

9080 

OUT 

820,0 

9090 

OUT 

820,18 

9100 

OUT 

820,0 

9110 

OUT 

820 , 21 

9120 

OUT 

820,0 

9130 

OUT 

820,8 

9140 

OUT 

820 , 0 

9150 

RETURN 

lOOOO 

1  REM 

SUB 

TO 

TURN  ON  9,10,16,18,22 

10040 

OUT 

821,255 

1007C 

OUT 

820 , 9 

loose 

OUT 

820,0 

1009C 

OUT 

820,10 

10100 

OUT 

820 , 0 

131 

10110  OUT  820,16 
10120  OUT  820,0 
10130  OUT  820,18 
10140  OUT  820,0 
10150  OUT  820, 22 
10160  OUT  820,0 


10161  REM  LINE  TO  KEEP  8  OFF 

10162  OUT  821,0 

10163  OUT  820,8 

10164  OUT  820,0 
10170  RETURN 

10500  REM  SUB  TO  TURN  OFF 

10510  OUT  821,0 

10511  REM  LINE  TO  KEEP  8  OFF 

10512  OUT  820,8 

10513  OUT  820,0 
10520  OUT  820,21 
10530  OUT  820,0 
10550  OUT  820,10 
10560  OUT  820,0 
10570  OUT  820,16 
10580  OUT  820,0 
10590  OUT  820,18 
10600  OUT  820,0 
10610  OUT  820,22 
10620  OUT  820,0 
10630  OUT  821,255 
10640  OUT  820,20 
10645  OUT  820,0 
10650  RETURN 

11000  REM  SUB  TO  TURN  OFF 

11010  OUT  821,0 

11020  OUT  820,18 

11030  OUT  820,0 

11040  OUT  820,10 

11050  OUT  820,0 

11060  OUT  820,16 

11070  OUT  820,0 

11080  OUT  820,22 

11090  OUT  320,0 

11100  OUT  821,255 

11110  OUT  820,20 

11120  OUT  820,0 

11130  OUT  820,11 

11140  OUT  820,0 

11150  RETURN 

12000  REM  SUB  TO  TURN  OFF 

12010  OUT  821,0 

12020  OUT  820,20 

12030  OUT  820,0 

12040  OUT  821,255 

12050  OUT  820,17 

12060  OUT  820,0 

12070  OUT  820,21 

12080  OUT  820,0 


21,10,16,18,22  AND  TURN  ON  20 


18,10,16,22,  AND  TURN  ON  20,11 


20  AND  TURN  ON  21,17 


rsi^^^^i**^**-***-***-**^-***-***-^^*^*  *^^**1  intf*!  i‘l  1*1  ri  i“i  ■*!  Cl  i^i  i*!  ■*!  Hi  a‘i  >‘i  ■“!  ■  i 


12090 
13000 
13010 
13020 
130  30 
13040 
13050 
13060 
13070 
13080 
13090 
13100 
13110 
13120 
13130 
13140 
13150 
13160 
13170 
13180 
13190 
14000 
14010 
14020 
14030 
14040 
14050 
14060 
14070 
14080 
14090 
14100 
14110 
14120 
14130 
14140 
14150 
14160 
14170 
1  5000 
15010 
15020 
15030 
15040 
15050 
15060 
15070 
15080 
15090 
15100 
15120 
15130 
15140 
15150 
15160 


RETURN 

REM 

OUT  821 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
REM 

OUT  821 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
RETURN 
REM 

OUT  821 


OUT  821 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  320 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
RETURN 
REM 

OUT  821 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 
OUT  820 


SUB  TO  TURN  ON  19,23 

,255 
,19 
,0 
,  23 

,0 

LINES  TO  KEEP  INTAKE  VALVES  OFF 


SUB  TO  TURN  OFF  10  AND  TURN  ON  11,23,17,19,21,9 

,0 

,  10 

,0 

,255 

,11 

,  0 

,  23 

,  0 

,  17 

,0 
,  19 
,  0 
,  21 
,0 
,  9 
,0 

SUB  TO  TURN  OFF  EVERYTHING 


PROGRAM  II:  PULSE  METHOD 


INE 

40 

50 

60 

70-80 

90-220 

230-240 

250-380 

390-1430 

1440-1580 


DESCRIPTION 
Get  pulse  time 
Get  pulse  gap 
Disable  the  Multiplexer 
Close  gripper 

Lift  the  object  and  put  it  on  the  upper 

platform 

Open  gripper 

Return  to  the  initial  position 

Subroutines  for  operating  gripper  and  fingers 
Timer 


136 


10  REM  File  Name;  REPORT. BAS 

20  REM  PULSE  METHOD:  TYPICAL  MANEUVER 

30  REM  N — Delay;  Nl--Pulse  Time;  N2--Pulse  Gap 

40  Nl=20 

50  N2=60 

60  OUT  820,0 

70  Ml=8  :  GOSUB  390 

80  N=400  :  GOSUB  1440 

90  K=1  :  Ml=18  :  GOSUB  440 

100  K=8 

110  Ml=10  :  M2=20  :  M3=18  :  M4=16 
120  GOSUB  920 
130  K=5 

140  Ml=10  :  M2=16  :  M3=18  :  M4=21 

150  GOSUB  920 
160  K=1 

170  Ml=ll  :  M2=16  :  M3=22 
180  GOSUB  720 
190  K=3 

200  Ml=10  ;  M2=17  ;  M3=22  :  M4=18 

210  GOSUB  920 

220  N=200  :  GOSUB  1440 

230  Ml=9  :  GOSUB  390 

240  N=400  :  GOSUB  1440 

250  K=3 

260  Ml=16  :  M2=18  :  M3=10 

270  GOSUB  720 

280  K=2  :  Ml=20 

290  GOSUB  440 

300  K=4  :  Ml=20  ;  M2=ll 

310  GOSUB  560 

320  K=6 

330  Ml=ll  :  M2=17  ;  M3=21  :  M4=23 
340  GOSUB  920 
350  K=6 

360  Ml=17  :  M2=19  :  M3=21  :  M4=ll 
3''0  GOSUB  920 
380  END 

390  REM  SUB  1;  GRIPPER  OPERATE 

400  OUT  821,255 

410  OUT  820, Ml 

420  OUT  820,0 

430  RETURN 

440  REM  SUB  2:  ONE  FINGER  OPERATE 

450  FOR  L=1  TO  K 

460  OUT  821,255 

470  OUT  820, Ml 

480  OUT  820,0 

490  GOSUB  1490 

500  OUT  821,0 

510  OUT  820, Ml 

520  OUT  820,0 

530  GOSUB  1540 

540  NEXT 

550  RETURN 


560  REM  SUB  3:  TWO  FINGERS  OPERATE 

570  FOR  L=1  TO  K 

580  OUT  821,255 

590  OUT  820, Ml 

600  OUT  820,0 

610  OUT  820, M2 

620  OUT  820,0 

630  GOSUB  1490 

640  OUT  821,0 

650  OUT  820, Ml 

660  OUT  820,0 

670  OUT  820, M2 

680  OUT  820,0 

690  GOSUB  1540 

700  NEXT 

710  RETURN 

720  REM  SUB  4:  THREE  FINGERS  OPERATE 

730  FOR  L=1  TO  K 

740  OUT  821,255 

750  OUT  820, Ml 

760  OUT  820,0 

770  OUT  8 20, M2 

780  OUT  820,0 

790  OUT  8 20, M3 

800  OUT  820,0 

810  GOSUB  1490 

820  OUT  821,0 

830  OUT  8 20, Ml 

840  OUT  820,0 

850  OUT  820, M2 

860  OUT  820,0 

870  OUT  8 20, M3 

880  OUT  820,0 

890  GOSUB  1540 

900  NEXT 

910  RETURN 

920  REM  SUB  5:  FOUR  FINGERS  OPERATE 


930  FOR  L=1  TO  K 
940  OUT  821,255 
950  OUT  8 20, Ml 
960  OUT  820,0 
970  OUT  820, M2 
980  OUT  820,0 
990  OUT  820, M3 
1000  OUT  820,0 
1010  OUT  820, M4 
1020  OUT  820,0 
1030  GOSUB  1490 
1040  OUT  821,0 
1050  OUT  820, Ml 
1060  OUT  820,0 
1070  OUT  820, M2 
1080  OUT  820,0 
1090  OUT  820, M3 
1100  OUT  820,0 


1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 


OUT  820, M4 
OUT  820,0 
GOSUB  1540 
NEXT 
RETURN 
REM  SUB  6: 
FOR  L=1  TO 
OUT  821, 25? 
OUT  820, Ml 
OUT  820,0 
OUT  820, M2 
OUT  820,0 
OUT  820, M3 
OUT  820,0 
OUT  820,M4 
OUT  820,0 
OUT  820, M5 
OUT  820,0 
GOSUB  1490 
OUT  821,0 
OUT  8 20, Ml 
OUT  820,0 
OUT  8 20, M2 
OUT  820,0 
OUT  8 20, M3 
OUT  820,0 
OUT  820,M4 
OUT  820,0 
OUT  820, M5 
OUT  820,0 
GOSUB  1540 


FIVE 

K 


FINGERS  OPERATE 


OUT  820, 
OUT  820, 
OUT  820, 
OUT  820, 
OUT  820, 
OUT  820, 
OUT  820, 
OUT  820, 
OUT  820, 
OUT  820, 
GOSUB  1! 
NEXT 
RETURN 
REM  SUB 
FOR  1=1 
J  =  I  +  1 
NEXT 
RETURN 
REM  SUB 
FOR  1=1 
J  =  I  +  1 
NEXT 
RETURN 
REM  SUB 
FOR  1=1 
J  =  I  +  1 
NEXT 
RETURN 


TIMER 

N 


8:  TIMER 
TO  N1 


9:  TIMER 
TO  N2 


s  s 


